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' I review investigations of the quark confinement mechanism that have been carried 

out in the framework of SU(N) lattice gauge theory. The special role of Zjy center 
symmetry is emphasized. 
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'. 1 Introduction 

> ■ 

^ ■ When the quark model of hadrons was first introduced by Gell-Mann and Zweig in 1964, 
■ an obvious question was "where are the quarks?". At the time, one could respond that the 
quark model was simply a useful scheme for classifying hadrons, and the notion of quarks 
as actual particles need not be taken seriously. But the absence of isolated quark states 
became a much more urgent issue with the successes of the quark-parton model, and the 
introduction, in 1972, of quantum chromodynamics as a fundamental theory of hadronic 
physics [1]. It was then necessary to suppose that, somehow or other, the dynamics of the 
gluon sector of QCD contrives to eliminate free quark states from the spectrum. 

Today almost no one seriously doubts that quantum chromodynamics confines quarks. 
Following many theoretical suggestions in the late 1970's about how quark confinement 
might come about, it was finally the computer simulations of QCD, initiated by Creutz [2] 
in 1980, that persuaded most skeptics. Quark confinement is now an old and familiar idea, 
routinely incorporated into the standard model and all its proposed extensions, and the 
focus of particle phenomenology shifted long ago to other issues. 

But familiarity is not the same thing as understanding. Despite efforts stretching over 
thirty years, there exists no derivation of quark confinement starting from first principles. 



1 



nor is there a totally convincing explanation of the effect. It is fair to say that no theory of 
quark confinement is generally accepted, and every proposal remains controversial. 

Nevertheless, there have been some very interesting developments over the last few years, 
coming both from string/M-theory and from lattice investigations. On the string/ M-theory 
side there is an emphasis, motivated by Maldacena's AdS/CFT conjecture [3], on discov- 
ering supergravity configurations which may provide a dual description of supersymmetric 
Yang- Mills theories. While the original AdS/CFT work on confinement was limited to cer- 
tain supersymmetric theories at strong couplings, it is hoped that these (or related) ideas 
about duality may eventually transcend those limitations. Lattice investigations, from the 
point of view of hadronic physics, have the advantage of being directly aimed at the the- 
ory of interest, namely QCD. On the lattice theory side, the prevailing view is that quark 
confinement is the work of some special class of gauge field configurations — candidates 
have included instantons, merons, abelian monopoles, and center vortices — which for some 
reason dominate the QCD vacuum on large distance scales. What is new in recent years is 
that algorithms have been invented which can locate these types of objects in thermalized 
lattices, generated by the lattice Monte Carlo technique. This is an important development, 
since it means that the underlying mechanism of quark confinement, via definite classes of 
field configurations, is open to numerical investigation. 

Careful lattice simulations have also confirmed, in recent years, properties of the static 
quark potential which were only speculations in the past. These properties have to do 
with the color group representation dependence of the string tension at intermediate and 
asymptotic distance scales, and the "fine structure" (i.e. string-like behavior) of the QCD 
fiux tube. Taken together, such features are very restrictive, and must be taken into account 
by any proposed scenario for quark confinement. 

This article reviews some of the progress towards understanding confinement that has 
been made over the last several years, in which the lattice formulation has played an impor- 
tant role. An underlying theme is the special significance of the center of the gauge group, 
both in constructing relevant order parameters, and in identifying the special class of field 
configurations which are responsible for the confining force. I will begin by focusing on what 
is actually known about this force, either from numerical experiments or from convincing 
theoretical arguments, and then concentrate on the two proposals which have received the 
most attention in recent years: confinement by center vortices, and confinement by 
abelian monopoles. Towards the end I will also touch on some aspects of confinement in 
Coulomb gauge, and in the large Nc hmit. I regret that I do not have space here to cover 
other interesting proposals, which may lie somewhat outside the mainstream. 

2 What is Confinement? 

First of all, what is quark confinement? 

The place to begin is with an experimental result: the apparent absence of free quarks 
in Nature. Free quark searches are basically searches for particles with fractional electric 
charge [4, 5], and the term "quark confinement" in this context is sometimes equated with 
the absence of free (hadronic) particles with electric charge ±|e and ±|e. But this exper- 
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imental fact, from a modern perspective, is not so very fundamental. It may not even be 
true. Suppose that Nature had supphed, in addition to the usual quarks, a massive scalar 
field in the 3 representation of color SU(3), having otherwise the quantum numbers of the 
vacuum. In that case there would exist bound states of a quark and a massive scalar, which 
together would have the flavor quantum numbers of the quark alone. If the scalar were not 
too massive, then fractionally charged particles would have turned up in particle detectors 
(and/or Millikan oil-drops) many years ago. It is conceivable (albeit unlikely) that such a 
scalar field really exists, but that its mass is on the order of hundreds of GeV or more. If 
so, particles with fractional electric charge await discovery at some future detector. 

Of course, despite the fractional electric charge, bound state systems of this sort hardly 
qualify as free quarks, and the discovery of such heavy objects would not greatly change 
prevailing theoretical ideas about non-perturbative QCD. So the term "quark confinement" 
must mean something more than just the absence of isolated objects with fractional electric 
charge. A popular definition is based on the fact that all the low-lying hadrons fit nicely 
into a scheme in which the constituent quarks combine in a color-singlet. There is also no 
evidence for the existence of isolated gluons, or any other particles in the spectrum, in a 
color non-singlet state. These facts suggest identifying quark confinement with the more 
general concept of color confinement, which means that 

There are no isolated particles in Nature with non-vanishing color charge 

i.e. all asymptotic particle states are color singlets. 

The problem with this definition of confinement is that it confuses confinement with 
color screening, and applies equally well to spontaneously broken gauge theories, where 
there is not supposed to be any such thing as confinement. The definition even implies the 
"confinement" of electric charge in a superconductor or a plasma. 

Imagine introducing some Higgs fields in the 3 representation into the QCD Lagrangian, 
with couplings chosen such that all gluons acquire a mass at tree level. In standard (but 
somewhat inaccurate [6]) terminology, this is characterized as complete spontaneous sym- 
metry breaking of the SU(3) gauge invariance. As in the electroweak theory, where there 
exist isolated leptons, W's, and Z bosons, in the broken SU(3) case the spectrum will con- 
tain isolated quarks and massive gluons. But does this really mean that there is no color 
confinement, according to the definition above? Note that whether the symmetry is broken 
or not, the non-abelian Gauss law is 

V • ^« = -c«^^4£;^ - ^^t^<l> (1) 

where c"*'^ are the structure constants, {t^} the SU(3) group generators, and (f) represents 
the Higgs and any other matter fields. The right-hand side of this equation happens to be 
the zero-th component of a conserved Noether current 

■a ^ _^al>c^,^pc, _ i-^t^<l> (2) 

and therefore the integral of Jq in a volume V can be identified as the charge in that volume. 
Just as in the abelian theory, the integral of the electric field at the surface of V measures 
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Bound antiquark screens 
quark electric field 
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Higgs condensate screens 
quark electric field 



Figure 1: In both QCD with dynamical quarks, and in a theory with complete spontaneous 
symmetry breaking, the field of a static color source Q in the fundamental representation is 
completely screened. 



the charge Q"" inside: 



Now suppose that there is a quark in volume V, far from the surface. Since the gluons are 
massive, the color electric field due to the quark falls off exponentially from the source, and 
the charge in volume V, as determined by the Gauss law, is essentially zero. So the isolated 
quark state, if viewed from a distance much greater than the inverse gluon mass, appears 
to be a color singlet. 

What is going on here is that the color charge of the quark source is completely cancelled 
out by contributions due to the Higgs fields and the gauge field surrounding the source. The 
very same effect is found in the abelian Higgs model, which is a relativistic generalization of 
the Landau-Ginzburg superconductor. Here too, the Higgs condensate rearranges itself to 
screen the electric charge of any source. Because of this screening effect there are no isolated 
particles, in the spectrum of a spontaneously broken gauge theory, which are charged with 
respect to a generator of the broken symmetry. 

For QCD, the situation is as depicted in Fig. 1. We consider a heavy quark Q sitting 
at the origin. In ordinary, "unbroken" QCD, the color electric field of the heavy quark 
is absorbed by a light antiquark, or any other set of quarks or colored-charged scalars 
forming a color singlet bound state with the heavy quark. In the spontaneously broken 
case, where the gauge group is completely broken, the charge of the heavy quark is shielded 
by the compensating charge of the gauge and condensate fields. In either case the total 
color charge in volume V is zero; there are no asymptotic particles in the spectrum with 
non-vanishing color charge. 

The similarity of the broken and unbroken gauge theories is not an accident. Fradkin 
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and Shenker [7], in 1979, considered a lattice model 

-Sfs = PgJ2J2[^ [U.i^W^ + nHix + 9)Ul{x)] + c.c.} 

+ PHY.12{^^i^Wx)cf>{x + il)+ c.c.} , 101 = 1 (4) 

X fj, 

which interpolates between the Higgs {PhjPg oo) and the confinement {i3h,Pg ^ 0) 
limits. They were able to show that for a Higgs field in the fundamental representation of 
the gauge group, the two coupling regions are continuously connected, rather than being 
separated everywhere in coupling space by a phase boundary. That is consistent with the 
fact that there are only color singlet asymptotic states in both regimes. 

The absence of a phase separation may seem paradoxical, if we imagine that the gauge 
symmetry is spontaneously broken at large (3h,(^g and restored at small Ph,Pg- There is 
no real paradox, of course, but unfortunately the phrase "spontaneous breaking of gauge 
symmetry" , although deeply embedded in the lexicon of modern particle physics, is a little 
misleading. In fact there is no such thing as the spontaneous breaking of a local gauge sym- 
metry in quantum field theory, according to a celebrated theorem by Elitzur [6]. However, 
after removing the redundant degrees of freedom by some choice of gauge, there typically 
remain some global symmetries, such as a global center symmetry, which may or may not be 
spontaneously broken. In the case that there are matter fields in the fundamental represen- 
tation of the gauge group, there is no residual center symmetry, and the Fradkin- Shenker 
result assures us that there is no symmetry breaking transition of any kind which would 
serve to separate a Higgs from a confining phase. The notion that confining and Higgs 
physics, in theories with fundamental matter fields, are separated by a symmetry breaking 
transition serves only as a "convenient fiction" [8] . 

The fiction is convenient because the spectrum of theories with "broken" gauge symmetry 
is qualitatively very different from the QCD spectrum. This qualitative difference is not 
captured very well by the notion of color confinement as it was defined above, because that 
property is found in both the Higgs and the "confining" coupling regions. What really 
distinguishes QCD from a Higgs theory with light quarks is the fact that meson states in 
QCD fall on linear, nearly parallel, Rcggc trajectories. This is a truly striking feature of the 
hadron spectrum, it is not found in bound-state systems with either Coulombic or Yukawa 
attractive forces, and it needs to be explained. 

2.1 Regge Trajectories and Color Fields 

When the spin J of mesons is plotted as a function of squared meson mass m^, it turns 
out that the resulting points can be sorted into groups which lie on straight lines, and that 
the slopes of these lines are nearly the same, as shown in Fig. 2. These lines are known as 
"linear Regge trajectories," and the particles associated with a given line all have the same 
fiavor quantum numbers. Similar linear trajectories are found for the baryons, out as far as 
J = 11/2. 

This remarkable feature of hadron phenomenology can be reproduced by a very simple 
model. Suppose that a meson consists of a straight line-like object with a constant energy 
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Figure 2: Reggc trajectories for the low-lying mesons (figure from Bali, ref. [9]). 



density a per unit length, having a nearly massless quark at one end of the line, and a 
nearly massless antiquark at the other. The quark and antiquark carry the flavor quantum 
numbers of the system, and move at nearly the speed of light. For a straight line of length 
L — 2R, whose ends rotate at the speed of hght, the energy of the system is 

^ adr 



m — E 



TiaR (5) 



while the angular momentum is 

J 



i 



^ arv{r)dr 

_ 2 pR ar^dr 
RJo JT^ 



^naR' (6) 

which means that this very simple model has caught the essential feature, namely, a linear 
relationship between and J. From the particle data, the slope of the Regge trajectories 



Comparing m and J, we find that 
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is approximately 

a' ^-^^ 0.9 (8) 
27r(7 ^ ^ 

implying an energy /unit length of the line between the quarks, which is known as the "string 

tension" , of magnitude 

a ^ 0.18 GeV^ ^ 0.9 GeV/fm (9) 

Of course, the actual Rcgge trajectories don't intercept the x-axis at = 0, and the 
slopes of the different trajectories are slightly different, as can be seen from Fig. 2. But the 
model can also be modified by allowing for finite quark masses. Note that since a crucial 
aspect of the model is that the quarks move at (nearly) the speed of light, the low-lying 
heavy quark states (charmonium, "toponium", etc.), composed of the c.t.b quarks, would 
not be expected to lie on linear Regge trajectories. Another way of making the model more 
realistic would be to allow for quantum fluctuations of the line-like object in directions 
transverse to the line. Those considerations lead to (and in fact inspired) the formidable 
subject of string theory [10]. 

QCD can be made agree with the simple phenomenological model if, for some reason, 
the electric field diverging from a quark is collimated into a fiux tube of fixed cross-sectional 
area. In that case the string tension is simply 



/ 



(fx± ^E^-E" (10) 



where the integration is in a plane between the quarks, perpendicular to the axis of the 
fiux tube. The problem is to explain why the electric field between a quark and antiquark 
pair should be collimated in this way, instead of spreading out into a dipole field, as in 
electrodynamics, or simply petering out, as in a spontaneously broken theory. 

In fact, as already emphasized, the color electric field of a quark or any other color 
charge source does peter out, eventually. If a heavy quark and antiquark were suddenly 
separated by a large distance (compared to usual hadronic scales), the collimated electric 
field between the quarks would not last for long. Instead the color electric fiux tube will 
decay into states of lesser energies by a process of "string breaking" (Fig. 3), which can 
be visualized as production of light quark- antiquark pairs in the middle of the flux tube, 
producing two or more meson states. The color fleld of each of the heavy quarks is flnally 
screened by a bound light quark, as indicated in Fig. 1. This process also accounts for the 
instability of excited particle states along Regge trajectories. 

Pair production, however, is suppressed if all quarks are very massive. Suppose the 
lightest quark has mass rUq. Then the energy of a flux tube state between nearly static 
quarks will be approximately while the mass of the pair-produced quarks associated 
with string-breaking will be at least 2mq. This means that the flux tube states will be 
stable against string breaking up to quark separations of approximately 

(11) 

a 

This brings us to the following statement of the problem we are interested in: 
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Figure 3: String breaking by quark-antiquark pair production. 

The Confinement Problem (I) 

Show that in the hmit that the masses of all quarks go to infinity, the work required to 
increase the quark separation in a quark-antiquark system by a distance L approaches aL 
asymptotically, where a is a constant. 

2.2 A First Encounter with the Center 

It seems that in order to define "quark confinement," we must effectively remove quarks as 
dynamical objects from the gauge theory. In the infinite mass limit, of course, quarks do not 
contribute to any virtual process, including string-breaking processes. But the exclusion of 
quarks doesn't mean that all matter fields must be removed, or taken to the infinite mass 
limit; we have only to exclude those matter fields which can give rise to string breaking. 
The criterion here is group-theoretical: If it is not possible for an individual quark and some 
number of matter field quanta to form a color singlet, then it is also not possible for the 
matter field to give rise to string-breaking. Suppose, for example, that QCD included some 
scalar fields in the color 8 representation. It is not possible for a particle or particles in the 
8 representation to combine with a quark in the 3 or antiquark in the 3 representation to 
form a color singlet, and therefore there can be no string breaking of the type shown in Fig. 
3. 

On the other hand, a set of Higgs field in the 8 representation of color SU(3) can still 
break the symmetry in such a way that all the gluons acquire a mass. In that case only a 
finite amount of work is required to separate two massive quarks by an arbitrary distance, 
even in the rrig oo limit. QCD with finite-mass matter fields in the 3 representation is 
therefore quite different from QCD with finite-mass matter fields in the 8 representation. 
In the former case, according to the Fradkin-Shenker result, there is no phase transition 
from the Higgs phase to a confining phase, and the work required to separate quarks by 
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an infinite distance always has a finite limit. In the latter situation, there can exist a true 
phase transition between a confining phase, and a non-confining Higgs phase. Which phase 
is actually realized is a dynamical issue, and will depend on the shape of the Higgs potential. 

Matter fields in color representations which cannot give rise to string-breaking, such 
as the color 8 representation of SU(3), are known as fields of zero N-ality. N-ality, as it 
is referred to in the physics literature, or the representation "class", as it is known in the 
mathematical literature, refers to the transformation properties of a Lie group representation 
with respect to gauge group center. We recall that the center of a group refers to that set 
of group elements which commute with all other elements of the group. For an SU(N) gauge 
group, the center elements consist of all g G SU{N) proportional to the N x N unit matrix, 
subject to the condition that det{g) — 1. This is the set of N SU(N)elements group elements 

~ 1 

1 



/27rm\ 



[n 



0,l,2,...,iV-l) 



(12) 



These center elements form a discrete abelian subgroup known as Zj^. 

Although there arc an infinite number of representations of SU(N), there arc only a finite 
number of representations of Z]\f, and every representation of SU{N) falls into one of N 
subsets, depending on the representation of the Z^v subgroup in the given representation. 
Each representation in a given subset has the same N-ality, which is an integer k defined as 
the number of boxes in the corresponding Young tableau, mod N. Transformation by Zn G 
Zjv, for each representation of N-ality k, corresponds to multiplication by a factor exp( ^^^" ). 
Group representations of N-ality k = are special, in that all center transformations are 
mapped to the identity. 

It is easy to see that a quark of non-zero N-ality can never form a color singlet by binding 
to one or more particles of zero N-ality. According to the usual rules of group theory, the 
possible irreducible color representations of the bound states are formed by combining the 
boxes in the Young tableaux of the constituents. If only the quark has non-zero N-ality, 
then the N-ality of the bound state is identical to that of the quark. 

If follows that matter fields in N-ahty=0 color representations cannot cause string break- 
ing, and it is therefore unnecessary to take their masses to infinity, in order to properly define 
quark confinement. We can therefore restate the quark confinement problem a little more 
generally as follows: 



The Confinement Problem (II) 



Consider an SU(N) gauge theory with matter fields in various representations of the 
gauge group, and take the hmit that the masses of the non-zero N-ality matter fields go to 
infinity. Show that in this limit there exists a confining phase, in which the work required 
to increase the separation of a non-zero N-ality particle-antiparticle pair by a distance L 
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approaches aL asymptotically, where a is a (representation-dependent) constant. 



The next task is to formulate order parameters which can distinguish the confining phase 
of an SU(N) gauge theory from other possible phases 



3 Signals of the Confinement Phase 
3.1 The Wilson Loop 

We begin with the lattice SU(N) gauge theory Lagrangian containing a single, very massive 
(Wilson) quark field in a color representation r of SU(N) 

S = /3i:(l-^ReTr[t/(p)]) 

+ E I K + 4aM^)V'(^) - ^ E V^(^)(« + 7,)Uj[Hx)^{x + fi) \ (13) 

X [ ^ n=±l J 

where p denotes plaquettes, = —7^, < o; < 1, and f/^^^ is the link variable in 
representation r, with t/_^(x) = t/^(a;— /i). A Wick rotation to imaginary time is understood. 
Suppose we create a quark- antiquark pair at time t — Q separated by a distance it! along, 
say, the x-axis, and let this system propagate for a time interval T. In the absence of gauge 
fixing, the expectation value of a color non-singlet state will average to zero, so it is necessary 
to include a product of link variables (a "Wilson line" ) between the quarks in order to form 
a gauge-invariant creation operator: 

g(i) = V^(o,t)r nVw(m,t)^(i??,t) (14) 

n=0 

where V is some 4x4 matrix, constructed from Dirac 7 matrices, acting on the spinor 
indices. Then, by the usual rules of quantum mechanics in imaginary time. 



{Q\T)Qm 



Z,„.(o|Q^I")("k-"-n"0(/"IQ|o) 

(0|e-f^^|0) 



where indicates the Euclidean vacuum expectation value, and A£^„ = En — Eq. Now 
integrate out the quark fields in the path integral. The leading contribution, at large ruq, is 
obtained by bringing down from the action a set of terms ip{a + 'y4)Uilj along the shortest 
lines joining the quark operators in Q and , and these shortest lines form the timelike 
sides of an X T rectangle. This contribution corresponds to the two quarks propagating 
from t = to t = T without changing their spatial positions, and in the rUq — > 00 limit all 
other quark contributions are negligible by comparison. We then have 



(gt(r)Q(O)) = ^lDUD^Di;Q^{T)Q{0)e- 



10 



C(m, + 4a)-2^^ / DU Xr[U{R,T)]e-'^ 

C(mg + 4a)-2^iy,(i?,T) (16) 



where \J (R, T) is the path-ordered product of hnks along the rectangular contour with 
opposite sides of lengths R separated by time T, Xrio) is the group character (trace) of 
group element g in representation r, and C is a constant arising from a trace over spinor 
indices. The Wilson loop Wr{R,T) is defined as the expectation value of Xr\U{R,T)\, and 
Su is the Wilson action of the pure gauge theory. It follows that 

E Icnre-^^"^ ~ C{m, + Aar'^Wr{R, T) (17) 

n 

In this relation, A£'„ refers to the energy, above the vacuum energy, of the n-th energy 
eigenstate having a non-vanishing overlap with the state created by Q. These can be un- 
derstood as the flux tube cigcnstates, and as T ^ oo, only the lowest energy contribution 
AEmin contributes. Subtracting the ln(mg + 4«) self-energy terms, which are independent 
of R and therefore irrelevant for our purposes, the /^-dependent part A£'^j„ is contained in 
the quantity 

Vr[R) = — hm log 



T- 



Wr{R,T) 



(18) 



which will be referred to from here on as the static quark potential. 

The confinement problem, then, is to show that Vr{R) has the asymptotic behavior 

Vr{R) ~ arR (19) 

at large R, for non-zero N-ality representations r. Note that in the nig — > cso limit, this is 
an observable of the pure gauge theory. Put another way, confinement is a property of the 
gauge theory vacuum in the absence of matter fields. 

In general, though, the confinement criterion for non-zero N-ality allows there to be finite 
mass matter fields in zero N-ahty representations, as discussed in section 2. In fact, gluons 
themselves are zero N-ality particles. This means that gluons can cause string-breaking for 
heavy quark sources in zero N-ality representations. This is an important point, which we 
will return to in the next section. It means that the confinement criterion can only refer 
to quarks in non-zero N-ality representations; the color field of a zero N-ality source will 
be screened in both the Higgs and confined phases, rather than coUimated into a fiux tube. 
Only the string tension cr^ of non-zero N-ality sources can serve as an order parameter for 
the confined phase, in which 

(7r 7^ for all color charge sources of non-zero N-ality (20) 
3.2 The Polyakov Line 

As already noted, a gauge theory with Higgs fields in the adjoint (or other zero N-ality) 
representation can have distinct confining and Higgs phases. Deconfining phase transitions 
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Figure 4: The global center transformation. Each of the indicated links in the t-direction, 
at t = is multiplied by a center element z. The lattice action is left unchanged by this 
operation. 



can also occur in a pure gauge theory at finite temperature. Any lattice Monte Carlo 
simulation is a simulation at finite temperature, with the inverse temperature equal to 
the extent Lt of the lattice in the time direction. Numerical studies of finite temperature 
phase transitions are carried out on x Lt lattices, with Lt < Ls, and deconfining phase 
transitions are found, at fixed gauge coupling, as Lt is reduced below some critical value. 

Different phases of a statistical system are often characterized by the broken or unbroken 
realization of some global symmetry (it is impossible to break a local gauge symmetry, as 
we know from Elitzur's theorem). In an SU(N) gauge theory, with only zero N-ality matter 
fields, there exists the following global Zjv symmetry transformation on a finite periodic 
lattice (Fig. 4): 

Uo{x,to) ^ zUo{x,to) , z e Zn , all f (21) 

with all other link variables unchanged. The action is obviously invariant under this trans- 
formation, with factors z and cancelling in every timelike plaquette at t — to (a property 
which is only true in general if 2; is a center element). For the same reason, any ordinary 
Wilson loop is invariant under this transformation. 

However, not every gauge-invariant observable is unchanged by the transformation (21). 
An important observable which is sensitive to this transformation is a Wilson line which 
winds once through the lattice in the periodic time direction 

Fix) = IV [C/o(f, l)C/o(f , 2)...C/o(f , Lt)] (22) 

and which transforms under (21) as 

P(f) ^ zP{x) (23) 

This observable, which is simply a Wilson loop with non-zero winding number in the time 
direction, is known as a Polyakov Loop. The global symmetry (21) can then be realized 
on the lattice in one of two ways: 

(pf")) — / ^ unbroken symmetry phase 

I non-zero broken Zn symmetry phase 
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The relation of the Polyakov loop VEV to confinement is quite direct. A Polyakov loop can 
be thought of as the world-line of a massive static quark at spatial position x, propagating 
only in the periodic time direction. Then 

(P(f)) = e-^^^* (25) 

where Fg is the free energy of the isolated quark. In the confinement phase, the free energy 
of an isolated quark is infinite, while it is finite in a non-confined phase. Therefore, on a 
lattice with finite time extension 

unbroken Z^- symmetry <^ confinement phase 

In other words, confinement can be identified as the phase in which global center symmetry 
is also a symmetry of the vacuum [11]. The Polyakov loop is a true order parameter: zero 
in one phase, non-zero in another, which associates the breaking of a global symmetry with 
the transition from one phase to another. 

This fact provides further insight into the Fradkin-Shenker result. A gauge theory with 
matter fields in the fundamental representation, such as the action in eq. (4), is not invariant 
under global center symmetry. Since the symmetry itself doesn't exist at the level of the 
Lagrangian, there can clearly be no phase transition between its broken and unbroken 
realization, and therefore no transition from the Higgs to a genuine confining phase. On 
the other hand, for matter fields in the adjoint (or any zero N-ality) representation, the 
Lagrangian does have global center symmetry, and a distinct confinement phase can exist. 

The transformation (21) can be generahzed: It is an example of a wider class of singulair 
gauge transformations. We consider gauge transformations on a periodic lattice 

Uo{x, t) ^ g{x, t)Uo{x, t)g\x, t + 1) (26) 

which are periodic in the time direction only up to a Zn transformation: 

g{x,Lt + l)^z*g{x,l) (27) 

which affects Polyakov loops in the same way as (21). The particular transformation (21) 
is generated by the singular transformation 

In the continuum gauge theory in a finite volume, the singular gauge transformation is 
again periodic only up to a center transformation 

g{x,Lt)^z*g{x,Q) (29) 

and the spatial gauge potentials Ak{x) transform in the usual way. The AQ^x^t) potential, 
however, transforms in the following way: At t 7^ 0, L^, 

M^{x, t) = g{x, t)AQg\x, t) - -g{x, t)dQg\x, t) (30) 

9s 
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as usual, where Qs is the gauge couphng. However, at i = 0, L^, we define 



A'oi.x, 0) = g{x, 0)Ao(x, Q)g\x, 0) - lim —g{x, e)dog\x, e) 
= Ao{x,Lt) 



= g{x,Lt)Ao{x,Lt)g^{x,Lt) -lim— g{x,Lt- €)dog^{x,Lt- e) (31) 

What this definition does is to drop the delta function at t = 0, Lt which would normally 
be present in (30) due to the discontinuity in the transformation g{x,t). That means that 
a "singular" gauge transformation is not really a gauge transformation, and this should not 
be a surprise. If singular gauge transformations were true gauge transformations, then all 
gauge-invariant observables, including Polyakov loops, would be unaffected by the trans- 
formation. The term "singular gauge transformation" is therefore slightly misleading, but 
because of common usage it will be retained here. 

3.3 The 't Hooft Operator 

In ordinary electrodynamics, the Wilson loop holonomy 



U{C) = exp 



ie (p dx^Ai^{x) 
c 



(32) 



of a space-like loop C is simply the exponential of the magnetic fiux through a surface 
bounded by the loop C. Although $s 7^ imphes that there must be non-zero field strength 
somewhere on any surface bounded by C, it is possible that for a large loop this field strength 
is localized far from loop C itself, and that the field strength in the neighborhood of the 
loop is actually zero. A familiar example is the vector potential in the exterior region of a 
solenoid, of radius i?, oriented along the z-asds 

^""-9 {r>R) (33) 



27rer 



where the field strength vanishes outside the coil at r > R. Although the field strength 
vanishes, we know that the vector potential in this region can still affect the motion of 
electrons. This is the well-known Bohm-Aharonov effect, which makes use of the fact that 
for a loop C winding around the exterior of the solenoid, we have U{C) ^ 1. 

The exterior solenoid field can be expressed as a singular gauge transformation of the 
classical A^ — vacuum state, with the transformation given by 



g{r, 9, z, t) = exp 



(34) 



This transformation is obviously aperiodic around the loop C 

g{r, 9 = 27r, z, t) = c'^'^gir, 9^0,z, t) (35) 
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with the aperiodicity associated with an element e~**^ G U{1) of the center of the gauge 
group (for an abehan gauge group such as U(l), the center of the group is the group itself, 
but we are anticipating the SU(N) case). The solenoid field in the exterior r > R region, at 
^ 7^ 0, 27r is given by 

A^{r > R, e, z, t) = -'-g{r, 6, z, t)d^g\r, 9, z, t) (36) 
while at 6* = 0, 27r we define 

A^{r,0,z,t) = - lim ^g{r,e,z,t)df,g\r,e,z,t) 
= A^{r,27r,z,t) 

= -Ji^-5'(^27r-e,^,t)a^5(^(r, 27r-e,2;,t) (37) 

Once again, this definition has the effect of dropping the delta function which would normally 
arise from the derivative of g{r, 9, z, t) along the hypersurface ^ = of discontinuity. As 
pointed out in the last section, this means that the singular gauge transformation is not a 
true gauge transformation (otherwise it could not possibly affect a gauge-invariant operator 
such as a Wilson loop). In the i? — > limit, the singular gauge transformation creates a line 
of infinite field strength along the z-axis. This is the R — limit of the interior field of the 
solenoid. 

The vector potential (33), resulting from the singular gauge transformation (34), is 
certainly not the unique form for the exterior field of a solenoid, and can be altered by any 
non-singular gauge transformation. What is essential is the aperiodicity (35) along some 
hypersurface of dimension D — 1 (in this case the hypersurface at ^ = 0). The hypersurface 
itself carries no action, and its position, apart from its boundary, is not a physical observable. 
The boundary of the hypersurface, however, is the solenoid, which carries non- vanishing field 
strength. In this example, in the R = limit, the boundary of the hypersurface is the z-axis, 
which carries an infinite field strength. Any loop C which winds n times around the z-axis, 
i.e. which has linking number n with the (infinite or periodic) z-axis, has a holonomy 
which is altered by the singular gauge transformation in this way: 

U{C) e±^"*«t/(C) (38) 

with the sign in the exponent depending on the orientation of the loop relative to the 
direction of the B-field. 

The linking of Wilson loops and vortices is a crucial concept, and deserves further elab- 
oration. In D=2 dimensions, a loop can be topologically linked to a point, in the sense 
that such a loop cannot be moved far from the point without actually crossing the point 
(Fig. 5). Obviously this doesn't make sense in D=3 dimensions, where a loop and point 
can be moved apart without crossing. Likewise, in D=3 dimensions, a closed loop can be 
topologically linked to another closed loop, but there is no such linking in D=4 dimensions 
where any two loops can be separated without crossing. In D=4 dimensions, closed loops 
can be topologically linked to surfaces. To visuahze this last statement, suppose that a 



15 




(a) (b) (c) 



Figure 5: A loop topologically linked to: a) a point, D — 2; h) another loop, = 3; 4) a 
surface, D — A. 



loop C, embedded in a 3-dimensional subspace, is topologically linked in D=4 dimensions 
to a closed surface S. In the 3D subspace, S appears as a closed curve C If S and C are 
topologically linked in D=4 dimensions, then C and C are topologically linked in the D=3 
subspace (otherwise C could be moved arbitrarily far away from S, without crossing S, in 
the D=3 subspace). The general statement, in any number of dimensions, is that a closed 
loop can be topologically linked to a hypersurface of co-dimension 2. In our example, the 
singular gauge transformation (35) actually creates field strength along the z-axis at every 
time t, i.e. a surface of field strength in the z-t plane, and a loop winding once around the 
z-axis at a fixed time t is linked topologically to this surface. 

Now we generalize to SU(N), and again consider transformations of the gauge field which 
can alter loop holonomies, without changing the action in the neighborhood of the loop. In 
D=4 dimensions, let V3 denote a compact, simply- connected Dirac 3-voIume, whose closed 
boundary is the surface S. Consider any loop C topologically linked to S, where the loop is 
parametrized by x{t), t G [0, 1], and x(l) = x{0) G V3. Let the gauge field along the curve 
C be transformed by a singular gauge transformation with the discontinuity 

g(x{l)) ^ zg(x{0)) , z E Zj, (39) 

which means that 

U{C) z*U{C) (40) 

and the transformed gauge potential A^{x) is defined in the usual way, except that a delta 
function arising from the discontinuity of g{x{T)) on the hypersurface V3 is dropped. As 
in the abelian case, the singular transformation creates a surface of infinite field strength 
along S in the continuum gauge theory, which is known as a thin center vortex. For 
an SU(N) gauge group there are — 1 possible vortices, corresponding to the number of 
elements in Zj^ different from the identity. As in the abelian case, it is possible to modify 
the transformation (39) near S, and smear out the thin vortex into a surface-like region of 
finite thickness, and finite field strength. This is a "thick" center vortex. 

A particular example of a singular gauge transformation on the lattice is the transfor- 
mation 

Uy{x,yQ,XA_) ^ zUy{x,yQ,XA_) for x > Xq (41) 
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Figure 6: Creation of a thin center vortex. The shaded plaquette, and all other x-y plaquettes 
at sites (a:o, yo,S±) form the center vortex. The stack of vortex plaquettes lie along a line in 
D — 3 dimensions, or a surface in D — A dimensions. 



with all other links unchanged, as indicated in Fig. 6. The (half-infinite) Dirac volume in 
this case is all sites with x > xq, y = yo- AH x — y plaquettes on the surface S, parallel to 
the z — t plane at x — xo, y — yo, are transformed by a center element 

U(P) zU{P) (42) 

The surface is a thin center vortex. 

Following 't Hooft, we now go to a Hamiltonian formulation, and consider an operator 
B{C) which creates a thin center vortex at a fixed time t on curve C (cf. ref. [12] for 
an explicit construction of this operator). This means that the gauge fields at time t are 
transformed by a singular gauge transformation satisfying (39) on any curve C at time t, 
parametrized by x{t), which has linking number one with loop C . Then 

U{C')B{C) = zB{C)U{C') , z^Zn (43) 

Using only this commutation relation, 't Hooft argued in ref. [13] that only area-law 
or perimeter- law falloff for the U{C), B{C) expectation values is possible, and that in the 
absence of massless excitations, the simultaneous behavior 

{U{C)) ~ e-'^^^^^ and (5(C)) ~ e'^^^^^ (44) 

is ruled out, where P{C) is the loop perimeter. From this it follows that if there are no 

massless excitations, then B{C) is an order parameter for confinement, in the sense that a 
perimeter-law falloff {B{C)) ~ exp[— 6P(C)] implies that the theory is in the confinement 
phase, while an area-law falloff for {B{C)) implies spontaneous breaking of at least part of 
the gauge symmetry. 
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3.4 The Vortex Free Energy 

On a finite lattice, it is impossible to create a single vortex sheet winding through the 
periodic lattice along, say, the z-t plane, because it requires a half-infinite Dirac volume. 
Instead, vortices which are closed by lattice periodicity can only be created or destroyed in 
pairs. 

There is however, a trick (due to 't Hooft [14]) known as "twisted boundary conditions," 
which forces the number of vortices winding through the periodic lattice to be odd, rather 
than even. Consider an SU(2) lattice gauge theory in D=4 dimensions, and imagine changing 
the sign P — > —(3 of the coupling on the set of x-y plaquettes at xq, yo and all z, t. It is not 
hard to see that the minimal action configuration of such a modified action, on a lattice 
which is infinite in the x-direction, is gauge equivalent to 

Uy{x,yo,x±) = —I for x > xq (45) 

with all other links equal to the identity matrix / {x± denotes axes perpendicular to the x, y 
directions). In this configuration the plaquettes at negative coupling are equal to — Tr [I], 
with all other plaquettes equal to -|-Tr [I] . The negative Uy links are located in a half- infinite 
Dirac volume. On a finite lattice, however, the Dirac volume cannot extend indefinitely in 
the positive x-direction, and must end on a center vortex sheet, which winds through the 
lattice in the z and t directions. There is no need for this vortex sheet to be thin, which costs 
a great deal of action at large /3. The action can be lowered if the vortex sheet has some 
finite thickness (just how thick is a dynamical issue at the quantum level, see section 6.6 
below). The upshot is that twisted boundary conditions, implemented by setting [3 —(3 
on a "co-closed" set of plaquettes (x-y plaquettes on a particular z-t plane closed by lattice 
periodicity), requires that there are an odd number (at least one) of thick center vortices in 
the periodic lattice. 

The magnetic free energy of a Z2 vortex in a finite volume V — L^LyLzLt is defined as 
the excess free energy of a periodic lattice with twisted boundary conditions, as compared 
to a lattice without twist, i.e. 

^-Fmg ^ _Z. (45) 
Z+ 

where Z± indicates the partition function with normal (-I-) and twisted (— ) boundary con- 
ditions. The "electric" field energy is defined by a Z2 Fourier transform 

-F^, sr^ Zz 

z=± ^+ 

= l-g-'^-s (47) 

Let C be a rectangular loop of area A{C). The following inequality was proven by Tomboulis 
and Yaffe [15]: 

A(C) 

{Tr [U{C)]) < {exp[-Fe,]} ^^^^^ (48) 

Therefore, if the free energy of a thick vortex surface running in the periodic z-t directions 
has an area-law falloff with respect to the cross-sectional x-y area 

Fmg = cLzLte-^'-^'^y (49) 
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Figure 7: The behavior of the vortex free energy vs. lattice extension (same in all directions). 
Prom Kovacs and Tomboulis, ref . [16] . 



then this is a sufficient condition for confinement, because it implies an area-law bound for 
Wilson loops. 

The ratio has been calculated numerically in SU(2) lattice gauge theory via 

lattice Monte Carlo, by Kovacs and Tomboulis in ref. [16]. Their result is shown in Fig. 7. 
The rapid rise of the ratio to unity, with increasing lattice extension in the x, y directions, 
is consistent with confining behavior (49) of the vortex free energy. The twist trick can also 
be used to insert center vortex flux through any plane, or several planes simultaneously [17]. 

The drop in vortex free energy with lattice size is naturally associated with the vortex 
"spreading out" to its natural thickness as lattice size increases. It is possible to insert 
constraints which prevent this spreading out [18], and then the behavior (49) is lost.^ Finite 
temperature studies involving the vortex free energy are also relevant [19, 20, 21] (see section 

6.4 below), because they involve "squeezing" vortices in one direction. 

3.5 Confinement Without Center Symmetry? 

All four of the confinement criteria listed above 

related constraint can be imposed on SU(2) Wilson loop operators, such that Tr[[/(C)] is effectively 
positive at weak couplings, and this very strong constraint on the loop removes the area law falloff [22]. That 
result is not unexpected, since the rapid falloflF of the Wilson loop is due to strong cancellations between 
positive and negative contributions to the loop expectation value. 
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1. area law for Wilson loops, 



2. vanishing Polyakov lines, 

3. perimeter-law behavior for 't Hooft loops, 

4. exponential area falloff of the vortex free energy, 

can only be satisfied if the Lagrangian is invariant under the global center symmetry of 
eq. (21). If, on the other hand, the theory contains matter fields of non-zero N-ality which 
completely break the center symmetry, then all Wilson loops have perimeter law falloff, all 
Polyakov loops are finite, and the Dirac volumes of center vortices carry an action density 
which is infinite in the continuum limit. We may ask whether it is possible to find some other 
observable or order parameter, not listed above, which would distinguish a "confinement" 
phase from a Higgs phase even when the global center symmetry is explicitly broken by 
matter fields in the fundamental representation of the gauge group. 

The Fradkin-Shenker result argues very strongly against this possibility. If there were 
some order parameter, call it Q, whose expectation value could distinguish between the 
confined phase and the Higgs phase of the theory, then there would have to be a fine 
of transitions completely isolating these regions from one another in the coupling-constant 
phase diagram. But there is no such line of transitions and, in consequence, no such operator 
Q, at least in a gauge theory with scalar fields in the fundamental representation of the 
gauge group. An order parameter distinguishing between the Higgs and Coulomb phases 
is a possibility, but an order parameter distinguishing between the Higgs and confinement 
phase appears to be ruled out. 

A case in point is the operator suggested by Fredenhagen and Marcu (FM) [23]. The 
idea behind their proposal is that if there exist dynamical matter fields, then one might 
redefine "confinement" to mean string breaking, rather than a linearly rising static potential. 
Consider the state 

I*..) - ; ] y.<t>\{x)u{c,y)Uy)\'^o) (50) 

^JW{R,R) i 

where the sum runs over all matter fields in the fundamental representation, C^y is the 
contour shown in Fig. 8, and \l'o is the vacuum. Sites x and y are at equal times, and 
we define R = \x — y\ as the spatial separation. The Wilson loop W{R,R) runs over a 
rectangular R x R contour, and the factor l/\/W is introduced so that the norm of \^xy) 
is of order one. 

The FM order parameter is defined as 



p = lim {^o\-^xyj 

it— ♦OO 

Ei(0l(x)C/(C,,)0,(y)) 



= lim 



R^oo W{R, R) 



(51) 



This order parameter is sensitive to the origin of the perimeter law behavior of Wilson loops. 
If the perimeter-law falloff of a Wilson loop is dominated by charge screening due to matter 
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Figure 8: Contour for the Fredenhagen-Marcu operator. 



fields, tfien the numerator and denominator in eq. (51) are comparable, and the ratio p is 
non-zero in the large R limit. This is the FM criterion for confinement. On the other hand, 
if the perimeter-law falloff of the Wilson loop is independent of charge screening, then the 
denominator may be much larger than the numerator, and p — in the limit. 

It should be clear that this is really a color screening criterion, rather than a confine- 
ment criterion. In the Fradkin-Shenker model, the FM criterion is certainly satisfied in the 
Pg:Ph <^ 1 regime, where the dynamics is "confinement-hke" , but it is also satisfied deep 
in the Higgs regime f3o, ^ 1, where screening also takes place. The FM criterion has, in 
fact, been applied numerically to the adjoint Higgs model (which actually has a transition 
between the confining and Higgs phases), with the matter field 4> at points x and y in the 
adjoint representation. It was found numerically in ref. [24] that the criterion p ^ was 
satisfied in both the Higgs and the confinement phases. 

The term "confinement" is sometimes taken to be synonymous with the absence of color- 
charged states in the spectrum, whether or not there exists a confining static potential. 
With this usage, a theory based on, e.g., the G(2) gauge group, which has a trivial center, a 
trivial first homotopy group (see below), and an asymptotically fiat static potential, belongs 
on the fist of confining theories [25]. However, terminology which essentially identifies 
confinement with color screening has drawbacks that have already been pointed out. It 
implies, for example, that electric charge is "confined" in an ordinary superconductor. In 
a gauge theory with scalars in the fundamental representation, it implies that there is 
confinement deep in the Higgs regime. In a gauge theory with scalars (and other matter 
fields) only in the adjoint representation, which really does have a transition between the 
Higgs and confinement phases, this use of the word "confinement" means that both phases 
are confining. In theories without a high-temperature deconfinement transition, the theory 
at high temperature would have to be taken as "confining" if the low temperature phase, 
which fulfils the FM criterion, is deemed to be in a confining phase. 

We conclude that while the FM operator is a good order parameter for color screening, 
there is a legitimate distinction to be made between the screening of charge by, e.g., a 
condensate of some kind, and confinement by a linear potential. This distinction means 
that confinement and color screening are not truly synonymous. The transition from a 
confining to a screened potential is always associated with the breaking of a global center 
symmetry, and in the absence of a non-trivial center symmetry, only screened or Coulombic 
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potentials can be realized. One would therefore expect that this symmetry is also important 
in understanding the origin of the confining force. 



3.5.1 The Case of SO (3) Gauge Theory 

There is still one puzzle: gluons actually belong to the adjoint representation of the gauge 
group, and it would seem to make no difference at all, in the continuum limit, whether 
the gauge group is SU(N) or SU{N)/Zj^. Both groups have the same Lie algebra, but in 
the latter case the center of the gauge group is trivial. Strictly speaking, the SU{N)/Zn 
theory is non-confining, simply because it is impossible to introduce color sources of non- 
zero N-ality. The static potential therefore tends to a constant at large distances. Still, as 
noted in ref . [26] , the presumed universality of the continuum fixed point suggests that the 
SU(N) and SU{N)/Zj^ theories should be in some sense equivalent in the continuum. For 
example, one might expect a finite temperature "deconfinement" transition in both cases, 
and whatever extended objects dominate the vacuum of the SU(N) gauge theory, in the 
continuum limit, should also dominate the vacuum of the SU{N)/Z]^ theory. 

This puzzle has recently been addressed by de Forcrand and Jahn in ref. [26], who find 
that the center vortex free energy remains a good order parameter for, e.g., the confinement- 
deconfinement transition both in SU(2) and in SO(3)=SU(2)/Z2 gauge theories. The key 
point is that center vortex creation is well defined in both SU(N) and SU{N)/Zn theories. 
Consider a path in the SU(N) group manifold, parametrized by r G [0, 1], which is traced by 
a singular gauge transformation g^r) around a closed loop x{t) in Euclidean space . On the 
SU(N) manifold, a vortex-creating transformation is discontinuous, i.e. g{l) — zg{0). The 
same transformation, mapped to the SU(N)/Zjv group, traces a continuous but topologically 
non-contractible path on group manifold. In this way, the vortex-creating transformations 
have an unambiguous topological signature both in SU(N) and Sl](N)/ Z^. Formally, the 
definition of center vortices rests on the existence of a non-trivial first homotopy group of 
the gauge group modulo its center. This homotopy group, 7ri{SU{N)/Z^) — Zjv, is the 
same for both the SU(N) and S'U(N)/Zj^ groups. 

In SU(2) lattice gauge theory, it is possible to specify the number of center vortices 
piercing a plane, mod 2, by imposing twisted boundary conditions as discussed above. De 
Forcrand and Jahn point out that S0(3) lattice gauge theory can be reformulated as an 
SU(2) lattice gauge theory in which all possible twisted boundary conditions are summed 
over. Their construction is motivated by the formulation of lattice SU(2) gauge theory in 
terms of S0(3) and Z2 variables, which was developed by the authors of refs. [22, 27, 28, 29]. 
Instead of being imposed boundary conditions, the "twists" z^i, = sgn(?7^,^) are observables 
in SO (3) lattice theory, and are extracted in the fiu plane from the quantity 

Vnu^ ^ n ^SnTrF[t/(Px,Mi^)] (52) 

where Px,iJ.v is a plaquette at point x in the fiv plane. With the help of this observable, it is 
possible to define and calculate the vortex (or "twist" ) free energies 

i^(^) = -log||f| (53) 
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in SO (3) gauge theory as well as SU(2) gauge theory. 



4 Properties of the Confining Force 

Every theory of confinement aims at explaining the linear rise of the static quark potential, 
which is suggested by the linearity of meson Regge trajectories. The linear potential, how- 
ever, is only one of a number of properties of the confining force that a satisfactory theory 
of confinement is obhgated to explain. A complete list would include at least the following: 

• Linearity of the Static Potential 

• Casimir Scaling 

• N-ality Dependence 

• String Behavior: Roughening and the Liischer term 

Each item on the list is supported by strong evidence from lattice Monte Carlo simulations; 
the last two items are also bolstered by fairly persuasive theoretical arguments. We will 
consider each item in turn. 



4.1 Linearity 

There is a theorem which can be proven in lattice gauge theory, which says that the force 
between a static quark and antiquark is everywhere attractive but cannot increase with 
distance; i.e. 

The proof of this statement, which holds in any number of spacetime dimensions, requires 
nothing more than the refiection positivity of the lattice action, and a clever application 
of Schwarz-type inequalities to certain Wilson loops used in potential calculations [30]. 
Therefore the linearity of the static quark potential is in fact the limiting behavior; the 
potential is constrained to be increasing but concave downward. If the quark potential grew 
at a rate faster than linear, it would violate the above bound on the second derivative. 

It has been clear for a long time, from lattice Monte Carlo simulations, that the static 
quark potential is asymptotically linear at long distances. In practice, it is useful to replace 
the Wilson line operator Q{t) in cq. (14) with an operator which has a larger overlap with 
the ground state of the QCD flux tube. A number of different operators have been used, 
but a typical method [31] ("APE smearing", see also refs. [32]) is to first compute for each 
spatial link 

[/;(x,t) = [/,(x,t) 

+cYl + ^' t)Ul{yi + t, t) + Ul{x - k, t)Ui{x - k, i)C/fe(x -k + t, t)]{55) 
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Figure 9: The static quark potential in SU(3) lattice gauge theory. The data points are 
compared to perturbative calculations with different renormalization schemes; the dotted 
line is a fit to the string-inspired potential a — tt/ (12i?) + aR. From Necco and Sommer, ref. 
[33]. 



and then project the "thick" link variables fZ- back to the SU(N) group manifold. This 
procedure can be iterated, with c a fixed parameter. If parameter c and the number of 
smearing iterations are chosen optimally, then the overlap of Q'{t), constructed from thick 
hnks, onto the flux-tube ground state is large. Defining W{R, T) to be timehke Wilson loops 
with thick spacelike links, the quantity 



V(i?,T) = -log 



W{R,T + l) 
W{R, T) 



(56) 



converges rapidly to its T oo limit, which is the static quark potential V{R). 

Figure (9), taken from ref. [33], displays some typical results for the static potential 
V{R) in SU(3) lattice gauge theory, obtained from the lattice Monte Carlo procedure. In 
this figure the Monte Carlo data is compared to perturbative calculations of the potential, 
with differing renormalization schemes. The linear growth of the potential beyond 0.7ro is 
evident from the figure. 

The X-axis of the figure shows quark separation in units of the Sommer scale Tq ~ 0.5 
fm, which is often used in studies of the static potential. The Sommer scale is defined as 
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the separation at which [34] 

F{ro)rl = 1.65 (57) 
where F(r) is the force between static quarks in the fundamental representation. 

4.2 Casimir Scaling 

"Casimir scaling" [35] refers to the fact that there is an intermediate range of distances where 
the string tension of static sources in color representation r is approximately proportional 
to the quadratic Casimir of the representation; i.e. 

C 

Or = y^CFF (58) 
Up 

where the subscript F refers to the fundamental representation. This behavior was first 
suggested in ref. [36]. The term '"Casimir scaling" was introduced much later, in ref. [35], 
where it was emphasized that this behavior poses a serious challenge to some prevaihng 
ideas about confinement (see also the earlier comments in ref. [37]). 

There are two theoretical arguments for Casimir scaling: dimensional reduction [38, 39] 
and factorization in the large-N limit [40]. We begin with factorization. The trace of a 
Wilson loop holonomy U{C), in a representation r of SU(N) gauge theory, is the group 
character Xr[U{C)], and this can always be expressed in terms of a sum of products of the 
group character in the defining representation F and its conjugate. For a large number N 
of colors, only the leading term 

Xr{g)-^XFi9)x*F\9)+OiN-') (59) 

is important. Here n+m <^ N is the smallest integer such that the irreducible representation 
r is obtained from the reduction of a product of n defining ("quark") representations, and 
m conjugate ( "antiquark" ) representations. Large-N factorization tells us that if A and B 
are two gauge-invariant operators, then in the — > oo limit 

{AB) = {A) {B) (60) 

Applying this property to Wilson loops in representation r, we have 

Wr{C) = {Xr[U{C)]) 

{Xf[U{C)xT[U{C)]) 

W]^+"' (61) 



From this it follows that at large N 



{n + 'm)aF (62) 



which is precisely Casimir scaling in this limit. Therefore Casimir scaling is exact, out to 
infinite quark separation, at N = oo. In particular, as first noted in ref. [37], the string 
tension of the adjoint representation must be twice the string tension of the fundamental 
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Figure 10: Numerical evidence for Casimir scaling in SU(3) lattice gauge theory. The solid 
lines are obtained from a fit of the potential in fundamental representation, multiplied by a 
ratio of quadratic Casimirs CrjCp. From Bali, ref. [42]. 



representation in the large N limit. Assuming that there is nothing singular or pathological 
about the large-N limit, it follows that some approximate Casimir scaling, at least up to 
some finite range of distances, should be observed at finite A^. In fact, from the strong- 
coupling expansion in lattice gauge theory, one finds for square L x L Wilson loops in the 
adjoint representation that 



where the strong-coupling diagram contributing to the perimeter-law falloff is shown in Fig. 
11. For sufficiently large L the second term dominates, and therefore the asymptotic string 
tension is zero. At iV ^ cxd only the first term is important, and the string tension in the 
adjoint representation is a a = ^ap, in agreement with Casimir scaling. At finite values of 



Wa[C] = N^e-^''^^' + e-^^""^^ 



(63) 
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Figure 11: Strong-coupling diagram contributing to the perimeter-law contribution to a 
Wilson loop in the adjoint representation. 



N, and strong-couplings, we would have aA ~ 2(7^ up to the distance 



L = W — log(Ar) + 16 + 4 

V C^F 



(64) 



which increases logarithmically with N at large N. 

The second argument for Casimir scaling is the (supposed) property of dimensional 
reduction at large distance scales. One argument for this property goes as follows [38] : The 
expectation value of a planar, spacelike Wilson loop in D — A dimensions can be expressed 
in terms of the vacuum wavefunctional 



with 



^0 = exp[-R[A]] 



(65) 



(66) 



where R[A] is some gauge- invariant functional of the spatial components Ak{x) of the gauge 
field at a fixed time. Suppose that R[A] can be expanded in a power series of the field 
strength and its covariant derivatives 



R[A] = J d^x[aTr [F^] + /3Tr [F^] + 7Tr [DFDF] + 



(67) 



For small amplitude, long-wavelength configurations, which are assumed to dominate the 
functional integral at long range, we then have 



"ifolA] ~ exp -a J (fxTr [F^] 



(68) 



Lattice Monte Carlo simulations have verified this form of the ground-state wavefunc- 
tional for two instances of Yang- Mills field configurations: non-abelian constant fields, and 
"abehan" ([Aj,74j] = 0) plane- wave configurations [41]. 
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Assuming the form (68) for the vacuum wavefunctional, we have that for sufficiently 
large Wilson loops 

= W^=\C) (69) 

and the calculation of the planar Wilson loop in a D=4 dimensional theory can be reduced 
to D=3. The reasoning can be repeated to express the expectation value of large planar 
Wilson loops in D = 3 dimensions to the corresponding D = 2 dimensional calculation, so 
we have 

W^/^=4(C) ~ iy/^=2(C) (70) 

But Wilson loops in D—2 dimensions can be calculated directly from perturbation theory 
(the answer is obtained, in an appropriate gauge, just by exponentiating one-gluon ex- 
change), and the resulting string tensions satisfy Casimir scaling exactly. On these grounds, 
one argues that string tensions in the D — 4 dimensional theory should also satisfy Casimir 
scaling. 

After the significance of Casimir scaling was re-emphasized in [35], numerical studies 
were undertaken by Bali [42] and by Deldar [43] for the SU(3) group, to check just how 
accurately this law is satisfied. It appears that string tensions satisfy Casimir scaling to 
quite a high degree of accuracy, as shown in Fig. 10, taken from ref. [42]. 



4.3 N-ality Dependence 

For finite N the Casimir scaling law must break down at some point, to be replaced by a 
dependence on the N-ality k of the representation 

ar = fik)aF (71) 

The reason for this N-ality dependence is color screening by gluons. Quarks in the adjoint 
representation, for example, have zero N-ality, and according to Casimir scaling 

2Ar2 ^ ^ 

The static quark potential for adjoint quarks separated by a large distance R would then 
be roughly V{R) = gaR- At some point it becomes energetically favorable for a pair of 
gluons to pop out of the vacuum and bind to each of the quarks, forming two "gluelumps" 
of mass Mgl- This is a string-breaking process, and if we only look at energetics it should 
be favorable at 

iJ.= ?^ (73) 

However, expressed in terms of either Feynman or lattice strong-coupling diagrams, string- 
breaking is a non-planar process, and its amplitude is suppressed by a factor of l/N^, as 
already seen in eq. (63) and Fig. 11. This explains why Casimir scaling is exact out to 
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infinite R in the N = oo limit, even if M^l is finite in the same hmit. But allowing for 
non-planar contributions at finite A^, we expect that Casimir scaling must break down for, 
e.g., square L x L Wilson loops at a length scale which grows only logarithmically with N. 

In general, consider a static quark- antiquark pair with the quark in representation r of 
N-ality k. As the separation R is increased, it can become energetically favorable to pop 
some number of gluons out of the vacuum to bind with the quark and antiquark charges. 
The energetically most favorable representation of the quark-gluon bound state will be the 
lowest dimensional representation r' with same N-ality as r, since gluons (themselves of 
zero N-ality) cannot change the N-ality of a source. But this means that the asymptotic 
string tension of the quark- antiquark pair is the same as that of a quark- antiquark pair in the 
lowest dimensional representations of the same N-ality, which implies that asymptotic string 
tension depends only on N-ality k, as indicated in eq. (71). This asymptotic string tension, 
depending only on N-ality, is known as the "k-string tension" a{k). The k-representations 
are the representations of lowest dimensionality of N-ality k. In terms of Young tableaux, 
these are denoted by a single column of k boxes. Since a charge in a /c representation 
cannot be screened to some lower representation by binding to gluons, it is reasonable to 
suppose that the fc-string tension is the same as the Casimir string tension, in which case, 
asymptotically, 

a. - a{k) = (74) 

for representation r of N-ality k. An alternative behavior is suggested by MQCD and softly 
broken M = 2 super Yang- Mills theory. 

Or ^ a{k) = —^ap (75) 
sm^ 

which is known as "Sine-Law scahng" [44, 45]. We note that for fixed k, Casimir and 
Sine- Law scaling are identical in the N ^ oo limit, 

a{k) = kap (76) 

Although we can be fairly sure that the asymptotic string tension depends only on 
N-ality, both from the color-screening argument and from explicit lattice strong-coupling 
calculations, it is actually rather difficult to see this dependence explicitly. The problem 
is clear from Fig. 10 above, where the color 8 (adjoint) representation of SU{S) appears 
to have a string tension a a which is 9/4 that of the fundamental representation out to the 
largest separation calculated. There is no hint of a crossover to 0"^ = 0, where the potential 
would go from linearly rising to flat. The reasons for this arc technical. The potentials in 
Fig. 10 were calculated at T -C i?, using the smearing technique discussed in section 4.1 
to increase the overlap to the lowest energy flux tube state. This method works very well 
for the fundamental representation. However, for the adjoint representation, it may be that 
the overlap of the "smeared" Wilson line with the two gluelump state is very small, and 
therefore one really has to go to time separations T ~ i? in order for this state to become 
dominant in the sum over states (15). As a practical matter, it is difficult to compute loops 
in the adjoint representation which are large enough to observe the string-breaking effect. 
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Figure 12: The adjoint and |x fundamental static potentials V{R) vs R, in D = 3 dimen- 
sional SU(2) lattice gauge at /5 = 6.0. The horizontal line at 2.06(1) represents twice the 
energy of a gluelump. From de Forcrand and Kratochvila, ref. [48]. 



This problem was originally overcome by operator-mixing methods, and crossover from 
Casimir-scaling to asymptotic behavior for adjoint representation charges was observed nu- 
merically in ref. [46]. In the last year, however, a clever algorithm invented by Liischer 
and Weisz [47] has greatly increased the attainable accuracy of Wilson loop calculations, 
and allows numerical calculation of much larger loops than was possible previously, even 
with no further increase in processor power. The algorithm was recently exploited by de 
Forcrand and Kratochvila [48] to compute V{R) from V{R, T) at up to T 2i?, and at 
the largest values of R the crossover behavior to = from the (approximately) Casimir 
scaling value is clearly seen. In Fig. 12 the adjoint representation potential is shown, for 
SU(2) lattice gauge theory at /3 = 6.0 in D = 3 dimensions, compared to the fundamental 
representation string tension multiplied by the Casimir ratio Ca/Cf = 8/3. The crossover 
from a linearly-rising, Casimir-scaling potential to the asymptotic fiat behavior is seen at a 
separation of 10 lattice spacings. Note also that, prior to color-screening, there is a small 
but significant deviation of the adjoint potential from exact Casimir scaling. 

Color- screening for representations beyond the adjoint have not yet been studied numer- 
ically. On the grounds of very general energetics arguments, however, the dependence of 
string tensions on N-ality alone, i.e. ar = cr{k) is sure to hold. 
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4.4 String Behavior 



Linear Rcgge trajectories are explained, as we have seen, by a model in which a meson is 
a rotating straight line, of constant energy density, running between a massless quark and 
antiquark. Allowing for quantum fluctuations of the line in the transverse direction leads 
to the Nambu-Goto action 

Sn = J— I rfV \ldet[g\ (77) 



where x^{a^^(j'^) are coordinates of the worldsheet swept out by the line running between 
the quark and antiquark, as it propagates through D dimensional spacetime. Wick rotation 
to Euclidean spacetime is understood, and Qab is the induced metric on the worldsheet 



da"" da'' 



(78) 



The Nambu-Goto action, although inspired from hadronic physics, is not an adequate 
fundamental theory of mesons. For one thing, quantization is only consistent in D=26 
dimensions, and for another, the lowest lying state is tachyonic. Nevertheless, if we view 
this action as an effective theory of the QGD flux tube (with an implicit high-momentum 
cutoff of some kind), then it is possible to explore the effect, on the static quark potential, of 
quantum fluctuations of the electric flux tube in the transverse directions. This investigation 
was first carried out by Liischer in ref. [49] (see also ref. [50]). 

The first step is to introduce static quarks by demanding that the string worldsheet is 
bounded by an i? x T loop, with T ^ R. Denote coordinates by — (xq, Xi, x±), and use 
the invariance of Nambu-Goto action under reparametrizations to set = xq, = xi, 
with the Rx T loop lying in the xq — xi plane at x± — 0. Then the static potential, in this 
model, is given by 

g-W ^ Jox^e-^^ 

= Jdx± exp £ dxol^''dxi^l + {dox^)^ + {dox^)^ + 0[{dx^y 



(79) 



Expanding the action to second order in the transverse fluctuations 

^-v{R)T ^ e-"^^ / Df^ exp -a T dx^ T dxi i]:{doX^f + hd^x^)- 
J Jo Jo \2 2 



-crRT 



det{-d^ - dl) 



{D-2)/2 



(80) 



where the determinant of the two-dimensional Laplacian operator, subject to Dirichlet 
boundary conditions on the R x T loop, is regulated and evaluated via standard zeta- 
function methods. The result (for T ^ R) is that 



det(-9^ - dl) 



{D-2)/2 



— exp 



D-2 nT 



12R 



(81) 
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which leads to the static potential 



^ aU - (82) 

The term proportional to 1/Ris known as the Liischer term, and its origin is completely 
different from the Coulombic term in 3+1 dimensional QCD due to one-gluon exchange at 
short distances. The Liischer term is due to quantum fluctuations of the string, rather than 
one-gluon exchange, and can viewed as the Casimir energy of the string arising from the 
Dirichlet boundary conditions. 

The bosonic string model also predicts how the thickness of a QCD flux tube depends 
on the separation R of the quarks. In QCD, the energy density due to the color electric field 
parallel to the flux tube is determined by the connected correlator of a timelike plaquette 
P parallel to the R x T loop C 

_ {Tr [U{C)]Tr [U{P)]) ,^....p... .... 

^ {Tr mC)]) ^^^^ 

and the thickness of the flux tube is determined from the falloff of £ in the directions x± 
transverse to the plane of the RxT loop. The analogous calculation, in string theory, is to 
calculate the loop-loop correlator G{Ci, C2, h) given by the functional integral over string 
worldsheets with two parallel loops Ci and C2, separated by a transverse distance h, as 
boundaries of the worldsheet. Loop Ci represents the large RxT loop, while C2 represents 
the probe plaquette. 

It is not essential for the loops to be rectangular, and it simplifies the calculation if both 
loops are circular and concentric. Providing the transverse displacement h is not too large 
compared to the radii Ri > R2 of loops Ci,C2, the leading contribution to the loop-loop 
correlation function can be computed semiclassically 

G{Ci, C2, h) = exp [-aAminiCi, C2, h)] (84) 

where Amin{Ci,C2,h) is the area of the minimal surface with boundaries at Ci,C2. The 
answer is found to be 



AminiCi, C2, h) ^ n{Rl - Rl) + --^—-h^ (85) 

\og[Ri/R2) 

Plugging this into (84), we see from the Gaussian falloff of G{Ci,C2,h) with h that the 
width of the flux tube grows logarithmically with the quark- ant iquark separation (which for 
this geometry is on the order of i?i). The semi classical calculation breaks down at some 
point, since at some h the surface degenerates into a sheet spanning the minimal area of 
Ci, a sheet spanning the minimal area of C2, with a tube of infinitesimal radius connecting 
to the two sheets. Beyond this point, a quantum treatment is required. One approach is to 
simply calculate, for the quantized Nambu string, the expectation value 

= {xl{xo = T/2, xi = R/2)) (86) 
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Figure 13: Coefficient of the Liischer term c(r) in D=3 and D—4 dimensions. The Nambu- 
Goto string action predicts c = —0.262 for D=4 and c = —0.131 for D=3. Prom Liischer 
and Weisz, ref . [54] . 



The answer can be expressed as a divergent sum over string modes. Cutting off the sum at 
some k = k^ax gives the result 

^ ^^\og{k^,,R) (87) 
(J 

As in the semiclassical calculation, the width of the flux tube is found to grow logarithmically 
with quark separation. 

Since the Nambu action is clearly not the correct description of quark-antiquark forces, 
because of the pathologies already mentioned, one might doubt that the Liischer term or the 
logarithmic growth of the flux tube width, derived from bosonic string theory, have anything 
to do with QCD. In recent years, however, developments in M-theory — speciflcally the 
AdS/CFT correspondence conjectured by Maldacena [3] — have lent strong support to the 
string picture of the QCD flux tube. The string description of the confining force, in the 
AdS/CFT approach, is not the old Nambu string. Instead it is a 10-dimcnsional superstring 
in a peculiar metric, which nonetheless describes the physics of D = 4 dimensional QCD. 
Both the Liischer term and the logarithmic growth of the flux tube are natural in this 
framework [51, 52]. The catch is that the AdS/CFT correspondence for QCD is (so far) 
restricted to strong-coupling calculations with a short-distance cutoff, much like strong- 
coupling lattice gauge theory. The AdS/CFT correspondence is beyond the scope of this 
review, and the reader is referred to ref. [3] for an exposition of the subject. 

One of the implications of a logarithmic growth in the width of the QCD flux tube with 
quark separation is that the width is infinite at R — oo. This is a clear prediction of string 



33 



models. In contrast, in strong-coupling lattice gauge theory the loop-plaquette correlation 
function (83) can be calculated to high orders in (3, and the width of the flux tube is 
found to be finite at i? = oo for small /3. However, the power series has a finite radius of 
convergence, and Pade and other extrapolation methods indicate that Woo — > oo at a finite 
(5 — Pc [53]. The phase transition at this point is known as the "roughening transition," 
which is known to occur in various models in statistical mechanics which involve the thermal 
fluctuations of a boundary. It therefore appears likely that = oo is a genuine feature of 
lattice QCD at weak couplings, in accordance with the prediction of string theory. 

Verification of the Liischer term via lattice Monte Carlo simulations is a difficult task, 
and until rather recently the numerical evidence favoring this term was unconvincing. It 
requires a very precise measurement of the static potential, at distances where the potential 
is almost linearly rising, to accurately detect the small deviation from linearity predicted by 
string theory. Recently, using a powerful algorithm for computing Polyakov line correlators 
[47] , Liischer and Weisz [54] were able to compute the static potential in SU(3) lattice gauge 
theory to very high accuracy, and extracted 

m = (88) 

from the second lattice derivative of V{R). The results for c{R) ys. R in D = 3 dimensions 
a,t P — 20, and D — A a,t P — 6.0, are shown in Fig. 13. In both cases, the asymptotic value 
of c{R) is consistent with the Casimir energy of bosonic string theory, i.e. the Liischer term 
in cq. (82). This result has been confirmed by the independent investigation of Juge et al. 
[55] (see also Necco and Sommcr [33] and Teper [56]). A common conclusion of all these 
investigations is that the Liischer term (which is evidence of string-like behavior) makes its 
appearance at a surprisingly short quark separation, on the order of 0.5 fm. 



5 The Center Vortex Confinement Mechanism 

In this and the following sections we explore some possible explanations of the confining 
force, beginning with the proposal which is most closely tied to the center symmetry 

[13, 57, 58, 59, 60]. 

There are several motivations for the center vortex mechanism. The first is the fact that 
the asymptotic string tension cr^, extracted from Wilson loops in the r representation of 
color SU(N), depends only on the N-ality kr of the representation; i.e. 

ar = a{kr) (89) 

In particular, ar — ii kr — 0. A second motivation is the fact that confinement is associated 
with the unbroken realization of a global center symmetry, as explained in section 3.2. 

Spontaneous or explicit breaking of this global symmetry results in a non-confining theory. 
Finally there is the idea, which is common to a number of proposals about the confinement 
mechanism, that vacuum configurations in pure SU(N) gauge theory can be decomposed into 
high-frequency, non-confining perturbations A^{x) around a smooth confining background 
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A^{x). The sub-leading, perimeter-law falloff is mainly due to fluctuations in ^4^, while the 
long- wavelength background is responsible for the area law.^ 

Let us focus in particular on large N-ahty = Wilson loops which have asymptotic string 
tensions 0"^ = 0. Somehow the long-wavelength confining fluctuations don't affect the 
zero N-ality loops, which implies 



Xr[U{C)] ^ Xr[U{C)] 



(90) 



where 



U(C) = Pexp 
U{C) = Pexp 
U{C) = Pexp 



c 



c 



dx^ A^,{x) 



dx'^ A^{x) 



(91) 



But how can eq. (90) be true, for any confining background The most straightforward 
possibility is that 

U{C) = Z{C)gU{C)g-^ (92) 

where Z{C) G Zjsi is a center element, and therefore U{C) = Z{C). The lattice configura- 
tions U^{x) which satisfy U{C) — Z[C) on any loop C can be identified uniquely: 



U^{x) = Zi_,{x)g{x)g ^{x + fl) , Zi_,{x) e Z 



N 



(93) 



These are the thin center vortices, created by singular gauge transformations. The field 
strength of such a configuration is singular on the vortex surface, but of course this singu- 
larity can be smoothed out, in a typical vacuum configuration, by the An{x) contribution. 

In the center vortex picture, confinement is simply a matter of short range correlations 
in the center flux on minimal surfaces. Suppose that we can subdivide any plane on the 
lattice into L x L squares of area such that if Ci and C2 are L x L loops around any 

two different squares, then 



{Z{C,)Z{C2)) ^ {Z{C,)) (Z(C2)) 



(94) 



Then consider a planar loop C of area E which encloses a number T^/T^min of the L x L 
squares Ej, each bounded by a loop Ci. Using the commutative property of the Z{Ci), and 
the assumed property of short range correlations among areas Ej 



ixMC)]) 



Xr 



n 

1=1 



z'^{a)iN 



xain] n 

i=l 

dr exp[—(7{kr)T] 



(95) 



^Note that is not necessarily entirely perturbative. It might include, e.g., short-range effects due to 
the finite thickness of monopoles, vortices, instantons, or any other solitons of importance. 
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where 



a{k) 



log (z'iQ) 



(96) 



is the string tension of the loop, depending only on the N-ality. 

5.1 Vortices as Local Minima of the Lattice Action 
5.1.1 A^>4 

We begin by showing that center vortices are local minima of the Wilson lattice action, and 
that they also must exist as local minima of long-range effective lattice actions. Confinement 
is then a matter of whether or not vortex surfaces percolate through the lattice volume. The 
proof for the Wilson action, which is due to Bachas and Dashen [61], is quite simple: Any 
thin center vortex in SU(A'") lattice gauge theory is gauge equivalent to 



with 



exp 



'27iinn{x) 



(97) 



N 



, n^(a;) = l,2,...,7V-l. 



Now write a small deformation of this configuration as 

U^{x) = Z^{x)V^{x), V^{x) = e*^-(^) 

with 

A^{x)^J2Al{x)ta, |^^(x)|«l. 

a 

Substituting (98) into the Wilson action 



2N 



J2(2N-TT[Up]-Tr[ul. 



we obtain 



•5 = ^ E (2A^ - ZpTr[Vp] - Z*pTr[V^]) 



2N 



(98) 
(99) 

(100) 
(101) 



Writing the product of l^-link variables around a plaquette in terms of a field strength 

1 



Vp = e'^^ ^In + iFp - ^F'p + 0{F'^) 



and also writing for Zp 



(2ninp\ 
Zp = exp [-^) 



the Wilson action becomes 



^ = 4 E (i - cos (^) ) . COS (^) nFl\ . am 



(102) 
(103) 

(104) 
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This action has a local minimum at Tr[Fp] = if the coefficient in front of Tr[Fp], at each 
plaquette P, is positive. That requirement is trivially satisfied at non-vortex plaquettes 
{np — 0). For vortex plaquettes (np 7^ 0), we have a positive coefficient, and hence a local 
minimum at Tr[F|] = 0, providing that 



cos 



f27inp\ . np 1 N-np 1 



This condition cannot be satisfied for A'" < 4. However, beginning with np = 1, 4 at A?" = 5, 
vortices appear as local minima of the Wilson action already at the classical level. 

This simple result can of course be extended [62] beyond the simple Wilson action, and 
therein lies its physical relevance. It is obvious that the action of thin vortices, stable or 
not, is singular in the continuum limit, and their contribution to the functional integral 
can be neglected. The configurations of physical interest are vortices which have some finite 
thickness d in physical units. In order to investigate the stability of vortices of some thickness 
d (or of lesser thickness, but including quantum fluctuations of the vortex surface up to that 
scale) , starting from a lattice action at spacing a, we imagine following the renormalization 
group approach, successively applying blocking transformations of the form 

^-s'lu'] ^Jdu S[U' - F(C/)]e-^[^l (106) 

where Ujj^ are links on the blocked lattice, and F{U) is the blocking function, until the lattice 
scale of the blocked action equals d. 

In this kind of approach, one usually assumes that only a few contours (plaquettes, 6 
and 8-link loops, etc.) are important in the effective action on the blocked lattice. Given 
an effective action at length scale d, it is possible to study "thin" (thickness= d) vortices at 
that scale. For simplicity, consider lattice actions which consist of plaquette (P) and 1x2 
rectangle (R) terms: 

Si^Co^(N- ReTr[U{P)]) + Ci (iV - ReTr[U{R)]) (107) 

P R 

This class of actions has been widely discussed in the lattice literature, and includes the 
tadpole- improved action [63], the Iwasaki action [64], and two-parameter approximations 
[65] to the Symanzik action [66] and the DBW2 action [67] . 

The thin vortex configuration, eq. (97), is easily shown to be a stable minimum of the 
action Sj providing both cq and ci are positive; the proof is as simple as in the case of the 
Wilson action. But for most improved actions of this type one has cq > and ci < 0, and 
the proof of stabihty is a httle more involved. It can be shown that the condition for vortex 
stability is, first, that the trivial vacuum U — I is a global minimum of Sj, which requires 

Co + 8ci > (108) 

and, second, that the vortex stability condition for the Wilson action, eq. (105) also holds. 
The proof of this statement can be found in ref. [62]. 

So it appears that the result first derived by Bachas and Dashen [61] is very robust, and 
holds for a large class of improved actions in coupling-constant space. It seems unlikely that 
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this result would be altered by adding a few more contours, as in certain proposed perfect 
actions [68]. If, in fact, center vortices are local minima of the effective action all along the 
renormalization trajectory, then their effects must become apparent at some scale. This is 
because the entropy factor increases with vortex surface area as 

exp[+const • (Vortex Area)] (109) 

while the Boltzmann suppression factor goes like 

exp[- (Vortex Area)/K^(ci)] (110) 

As d increases, so does ^^{d). Eventually entropy wins over action, and vortices at that 
scale will percolate through the lattice. 



5.1.2 N<4: 

We next consider N < A. Vortices are certainly not minima of the Wilson or two-parameter 
actions in this case, but there are good reasons to believe that they are local minima of the 
effective action at large scales. Suppose we have carried out RG blocking transformations 
up to some scale La which is well beyond the adjoint string-breaking scale. Define the 
probability distribution for loop holonomies U'{C) around contours C on the blocked lattice 

Pc{g) = {5[g.U\C)\) 

\ k Tk I 

= EEw^^.(c)x.J^1 (111) 

k rk 

where the sums are over N-ality A;, and over group representations of a given N-ality, and 
the Wr,.{C) are Wilson loop expectation values on the blocked lattice. The leading terms 
(largest Wr^{C)) come from the lowest k,N — k values, and, for each k, the representation 
of lowest dimension. Keeping only the leading terms atA; = 0,l,iV — Iwe have 

Pc{g) ^ l + WA{C)xA{g) + WF{C){xF{g)+c.c.) 

= 1 + e-'^P"™(^)-'^°XA(5) + e-'^^^'^^^(^)-'^^P^"'"(^)-'^i (xrig) + cc.) (112) 

where A, F refer to the adjoint and fundamental representations. 

Pc{g) is the probability density that a given loop holonomy has the SU(N) group value 
g. Because of color screening, the coefficient of XA{,g) is far greater than that of XF{,g) at 
large blocking scales. Since Xa(1) — Xa{z)-i this means that the loop probability density has 
local maxima at center elements. The only configurations which locally maximize all loop 
probabilities and are consistent with the Bianchi constraint are the center vortex configura- 
tions. We conclude that the effective action at this blocking scale has local minima at the 
vortex configurations. 
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This conclusion can be derived explicitly in strong-coupling lattice gauge theory, where 
the blocking transformation 

^ fj^ij Y[S[U'i, - (C/C/C/...C/)i,]e-^[^l (113) 
I' 

from a lattice of spacing a to a lattice of spacing La can be carried out analytically [69] . If 
g is the lattice coupling constant at lattice spacing a, one finds 

-|-larger adjoint loops with coefficients (114) 

The adjoint loops are non-planar contributions, and arise from the "tube" diagrams shown 
in Fig. 11. The non-local effective action can be expressed (with some additional complica- 
tions) as a local action containing a number of adjoint Higgs fields [69], but for simplicity 
we will simply consider the above action truncated to the leading plaquette terms, i.e. 

-S'trunAU'] = E{coReTV^[C/'(p')] + c,TrA[U'{p')]} (115) 

with 

Co ~ exp[— (7 Area(p')] , Ci ~ exp[—4a Perimeter (p')] (116) 
Consider small fluctuations around a vortex configuration (98), (102), (103) 

S'trunc = const + ^ E {co COS (^) Tv^F;,] + c^Tr^F^] } (117) 

It is clear that at large blocking L, Ci is much greater than Cq. Then S'trunc) the scale 
where Ci 3> Cq, is minimized at Fp/ = 0. This means that at this scale, vortex configurations 
are local minima of the coarse-grained action even for = 2, 3, 4. 

5.1.3 Center Monopoles 

Thin center vortices on a classical SU(N) vacuum background, which can be written in the 
form (97), can obviously be identified with configurations in Zj^ lattice gauge theory. In 
D dimensions these vortices form D — 2 dimensional surfaces. The Bianchi identities of 
Zi^ lattice gauge theory require that these surfaces are either closed, or else three or more 
surfaces may meet at D — 3 dimensional vertices, which are the monopoles of Zfq lattice 
gauge theory [70, 13]. 

It is easiest to visualize matters in D=3 dimensions. A Zj^ vortex pierces a given pla- 
quette p if Z{p) 7^ 1. We introduce a dual lattice whose sites lie in the center of unit cubes 
of the original lattice. Links of the dual lattice run through the middle of plaquettes on the 
original lattice. The vortex which pierces a set of plaquettes on the original lattice can be 
visualized as a path on the dual lattice. 
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Figure 14: Three vortices (solid lines) diverging from a monopole (solid circle) on the dual 
lattice, in Z3 lattice gauge theory. The vortices pierce plaquettes on the original lattice 
(dashed contours); oriented loops around these plaquettes, with normals in the direction of 
the vortices, have a value of exp[27ri/3]. 



Now consider a set of plaquette loops {pi,P2, --jPe} on ^ unit cube of the lattice, with 
the orientation of the loops chosen so that the normal to the plaquette points out of the 
cube. It is easy to verify the Bianchi identity 



In, e.g., Z3 lattice gauge theory, this identity allows three vortices to diverge from the dual 
lattice site at the center of the cube, as indicated in Fig. 14. The dual lattice site is then 
the position of a monopole. Of course this generalizes to higher A^: several vortices can 
diverge from a monopole providing the product of factors exp [27rm/A^] , associated with 
each vortex line, is equal to one. In D=4 dimensions vortices are surfaces on the dual lattice, 
joining at monopole worldlines. 

The local minima of the effective SU(N) lattice actions discussed above are simply the 
vortices and monopoles of a Zj\! lattice configuration. It has been argued recently that 
center monopoles play an important role in determining the precise N-ality dependence of 
asymptotic string tensions at large N [62]. 

5.2 Vortices as Local Minima of the One-Loop Effective Action 

In continuum SU(N) gauge theory, the only non-trivial classical solutions of the theory 
in D=4 Euclidean space are the instantons. The situation changes however, as shown 
by Diakonov and Maul [71], when quantum corrections at one-loop order are included. 
Diakonov and Maul calculate the Yang-Mills effective action, to one loop order, of the 
following configuration (in cylindrical coordinates) 



6 





i=l 




r 



(119) 
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where fi{r) is a profile of the vortex cross-section. Any large Wilson loop in the fundamental 
representation, located far from r = but winding n-times around the vortex, has the value 
^Tt[U{C)] = (—1)" in this configuration. This configuration is therefore a thick center 
vortex on a classical vacuum background. By the symmetry of the vortex configuration, the 
effective vortex action is proportional the extent of spacetime L^~'^ in the directions parallel 
to the vortex, and orthogonal to the (r, 0)-plane; i.e. 

3 vortex — ^-^f ^ (l^O) 

where S is the action density of the vortex in the transverse directions. Let tq be the scale of 
the vortex cross-section (e.g. /^(r) = exp[— ro/r]). Including quantum corrections up to one 
loop, Diakonov and Maul find a transverse action density for the vortex in D—A directions 

£ = ~^^-^,log[AQCDroB] (121) 

where A, B are dimensionless constants which depend on the precise functional form of //(r). 
Classically the vortex action is always minimized by sending the vortex thickness tq ^ oo. 
At the one loop level, however, the one-loop vortex action is minimized at a finite vortex 
thickness 

ro = (122) 



and action density at the minimum 



-^Ico^— (123) 



(where e = 2.718...). A similar analysis also finds a minimum thickness in = 3 dimensions, 
when one loop quantum corrections are included. The issue of vortex stability at one loop 
level is still not completely settled however; c.f. Bordag [72] for a discussion of this point. 



5.3 Locating Vortices in Monte Carlo Simulations 

It is a simple matter to create center vortices on a given background; this is accomplished by 
singular gauge transformations. The more difficult problem is to locate vortices in a given 
lattice configuration, but this can be managed by the technique of center projection in an 
adjoint gauge, which we now discuss. 

An adjoint gauge is a complete gauge-fixing condition for link variables Ua^i{x) in the 
adjoint representation, leaving a residual Zj^ symmetry for the links in the fundamental 
representation. The gauge-fixing condition may be that some functional of the links 7^[^7a] 
is zero or extremal, such as maximal center gauge, or it may fix the lattice to a unitary 
gauge of some kind, such as Laplacian center gauge. We can express the fundamental 
representation link variables in adjoint gauge as 

Ul^{x) = Z,{x)V,{x) (124) 
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where Z^{x) is the center element which is closest, on the SU(N) group manifold, to U^^{x). 
The mapping U^^ — > is known as Center Projection [73]. 

A center vortex is created, in an arbitrary configuration U, by making a singular gauge 
transformation. Call the result U'. The corresponding link variables in the adjoint repre- 
sentation, and U'^, are gauge equivalent. It follows that in an adjoint gauge, Ua and 
U'a are transformed into the same configuration 1/^ ■ The original fundamental link con- 
figurations U and U' are therefore transformed into configurations which can only differ by 
center elements, i.e. 

u';^{x) = z;(x)f/;^(x) 

= Z'^{x)Z,{x)V,{x) (125) 

where Z'^{x) is gauge-equivalent to the thin vortex created by the singular gauge transfor- 
mation. Therefore, each time a thin vortex is created in a given SU(N) configuration, the 
effect is to add a thin vortex to the corresponding center-projected configuration. For this 
reason, it is natural to interpret center-projected configurations Zfj_{x), of an SU(N) con- 
figuration in an adjoint gauge, as locating the position of center vortices in the full SU(N) 
configurations. The excitations of any Zn configuration are center vortices; the vortices 
of projected configurations are known as P-vortices, and the negative plaquettes in the 
projected configurations are known as P-plaquettes [73] . 

The weakness of this reasoning, which says that center projection in any adjoint gauge 
will locate the same vortices, is that it ignores the fact that center vortices in SU(N) gauge 
theories would have finite thickness, and different gauges may be more or less sensitive 
to this. In addition, many gauges suffer from the Gribov copy problem, and the center 
projection of different gauge copies may differ. Fortunately there are ways of checking 
whether P-vortices in a particular gauge actually locate thick vortices in the unprojected 
theory, by calculating the correlation of P-vortices with gauge-invariant observables. These 
correlations will be discussed in the next section. Some adjoint gauges which have been 
found useful for vortex finding in SU(2) lattice gauge theory are listed below. 

5.3.1 Direct McLximal Center Gauge [74] 

Direct maximal center gauge, which is simply the lattice Landau gauge in the adjoint repre- 
sentation, can also be regarded [75] as a best fit of a lattice configuration [/ by a thin vortex 
configuration U^{x) — g{x)Z^{x)g'^{x -\- /i). 

A typical thermalized lattice Un{x) generated in Monte Carlo simulations, if printed out, 
looks like a set of random numbers. But of course it is not random, and locally, at couplings 
/? 3> 1, the configurations are small fiuctuations around classical vacua 

U^ix)^gix)g''ix + il) (126) 

Suppose we ask for the best fit, to a given lattice configuration U^{x)^ by a pure gauge 
g{x)g\x + fi). We will consider, for simplicity, the SU(2) gauge group (the approach is 
readily generahzed to SU(N)). A best fit will minimize the distance on the SU(2) group 
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manifold 



4 = E Tr [{Uf^ix) - g{x)g\x + fi)) x (h.c. 



= J22Tr\l-g\x)U^{x)g{x + il) 



Define 



<^U^{x) = g\x)U^{x)g{x + i:i) 



Minimizing is the same as maximizing 



Y.n'u,{x)] 



(127) 



(128) 



which is the lattice Landau gauge. Fixing to lattice Landau gauge is therefore equivalent to 
the problem of finding the best fit of the given lattice to a vacuum configuration. 

In the same spirit, we may ask for the best fit of the lattice by a thin vortex configuration 



Wm(^) = 9i.x)Z^{x)g\x + /i) 



(129) 



In maximal center gauge, this is accomplished in two steps. First, since the adjoint repre- 
sentation is blind to Z^, find the best fit to the adjoint representation lattice by a (classical) 
vacuum configuration. The best fit minimizes 



4 = Y,'^A[{u^{x)-g{x)g'^{x + fi))xih.c. 
= Y,2TtA\l-g\x)U^{x)g{x + il) 



X,fJ, 



which is equivalent to maximizing the quantity 

R = J2TrA['U,{x)] 



(130) 



:i3i) 



X,IJ, 



The condition that it! is a maximum is the adjoint representation version of Landau gauge, 
and is also known as Direct Mctximal Center Gauge (DMC). 

The second step, for fixed g{x) determined by DMC gauge, is to choose Zfj^^x) so as to 
minimize the distance function in the fundamental representation 



^Y.T^t[{u^{x) - g{x)Z^{x)g\x + fi)) x (h.c.) 



and this is achieved by setting 

Z^{x) = signIVpC/^(x)] 
which is center projection in SU(2) gauge theory. 



(132) 



(133) 
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The original and gauge-transformed lattices can be expressed, respectively 

<^U^{x) = Z^{x)e''^^^^'^ (134) 

We interpret Z^{x) as the vortex background, and SA^^x) is the fluctuation around the 
background. DMC gauge finds the optimal Z^^x) minimizing 5An[x). 

Although this method of vortex-finding by a best fit procedure seems very natural, there 
are two shortcomings. First, in practice there is a Gribov copy problem. In general there is 
no technique for finding the global maximum of i?; instead there are various methods (over- 
relaxation, simulated annealing) which find local maxima. Secondly, one must be aware 
that the action density of a thin vortex is singular at P-plaquettes and therefore the fit to 
a thermalized lattice is very bad at those locations. In general 

^Tr[C/(p)] /3 > 1 (135) 

then 

^Tr[?7(p)] = Z{p)Ti[ W e*^^'*^^)] ^ 1 (136) 

links&p 

which implies 

-Tr[ [] e^^^''^^)] ^ -1 when Z{p) = -1 (137) 

links&p 

It follows that 8A^{x) is not small at these P-plaquette locations (and is in fact singular in 
the continuum limit). 



5.3.2 Indirect Maximal Center Gauge [73] 

Maximal abelian gauge (MAG) was introduced in ref. [76], in order to identify and study 
abelian monopoles in lattice configurations. In SU(2) lattice gauge theory, the gauge is 
defined as maximizing 

R^Y.'^[U^{x)a^Ul{x)a^] (138) 

which leaves a residual U(l) gauge symmetry. SU(2) link variables in this gauge can be 
decomposed as 

U^{x) = C^{x)D^{x) (139) 
where is the diagonal part of the link variable, rescaled to restore unitarity 



D, 



1 



i9, 



+ 



mJ22 



[U^]ii 

[U^]22 





e-*^" 



(140) 
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It can be seen that Dn[x) transforms like a U(l) gauge field under the residual U(l) gauge 
symmetry. 

Indirect maximal center (IMC) gauge begins with maximal abelian gauge fixing, and 
then uses the residual U(l) gauge symmetry to maximize 

R = Y.cos\9^{x)) (141) 

leaving a residual Z2 gauge symmetry. Vortex configurations are then identified by projec- 
tion 

Z^(a;) = sign [cos (^^(x))] (142) 

and this procedure can be interpreted as a best fit of vortex configurations to the abelian 
configuration D^{x) extracted from Un{x). 

Indirect center gauge has the advantage that both abelian monopoles and vortices can 
be identified and correlated. It has the same disadvantages as maximal center gauge; i.e. 
there is a Gribov copy problem, and the best fit is a bad fit at P-plaquette locations. 

5.3.3 Laplacian Center Gauge [77] 

Consider a Yang-Mills theory with two Higgs fields 0i, 02 in the adjoint representation. A 
unitary gauge, fixing color components 0i(a;) = p{x)5^^, leaves a residual U(l) symmetry. 
This residual symmetry can be used to set color components 02 (^) — 0; 02 (^) > leaving a 
remnant Z2 symmetry. The idea of Laplacian Center gauge, in pure Yang-Mills theory with 
no elementary scalars, is to replace the scalars 0i,2 with the two lowest-lying eigenstates (no 
sum over a, which labels eigenstates) 

^A,,(x,t/)/;(t/) = AX(^) (143) 
y 

of the lattice Laplacian operator in adjoint representation 

Aij(a;, y) = - ^ {[UAn{x)]ij5y,:c+ij. + [Uau.{x - fi)]jidy^x-f, - 25xy5ij) (144) 

There are standard numerical algorithms for obtaining these eigenstates. This construction 
has the great advantage that there is no Gribov copy problem. 

The use of eigenstates of the Laplacian operator for gauge fixing was originally proposed 
by Vink and Wiese [78] , in order to define a version of Landau gauge which avoids the Gribov 
copy problem. The use of the single lowest-lying eigenstate of the Laplacian operator in 
adjoint representation, to reduce the gauge symmetry from SU(2) to U(l), was proposed 
by van der Sijs [79], and is known as Laplacian abelian gauge, de Forcrand and co-workers 
[77] took this idea one step further, and used a second eigenstate of the adjoint Laplacian 
to also fix the U(l) symmetry, leaving only a remnant center symmetry. The extension of 
Laplacian center gauge to SU(N) theories is straightforward, and is described in ref. [77]. 

Langfeld et al. [80] have employed a combination of the LC and DMC gauges, first fixing 
to a unique configuration using the Laplacian center gauge, and subsequently relaxing to 
a nearby local maximum of the direct maximal center gauge. This amounts to selecting a 
particular Gribov copy of the DMC gauge, and is similar to the next center gauge we shall 
discuss. 
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5.3.4 Direct Laplacian Center Gauge [81] 

The direct laplacian center (DLC) gauge begins from a Gribov-copy free modification of 
adjoint Landau gauge. The idea is to first maximize 

Rm = J2T^t[M'^{x)Ua^{x)M{x + jl)] (145) 

with UAfj, in the adjoint representation. However, instead of imposing the condition (as 
in DMC gauge) that M{x) is an orthogonal matrix, we require only that M{x) is a real 
symmetric matrix with the orthogonality condition only imposed on average: 

i$:M^(x)M(x) = 7 (146) 

where V is the lattice volume. 

This problem is solved by finding the three lowest eigenvalues, and corresponding eigen- 
functions, of the lattice Laplacian eigenvalue equation (143). Then 

Mai,{x) = f,{x) (147) 

is the matrix field maximizing Rm, where fl,fi,ff are the three lowest-lying eigenstatcs. 
However, Mab is not an SO (3) gauge transformation, so we next map Mab{x) to the nearest 
SO (3) matrix- valued field 

[9A{x)]i, = fiix) (148) 
which also satisfies a Laplacian equation 

Y,^ijM~f!{y) = Kc{x)n{x) (149) 
y 

This is done by polar decomposition 

M{x) = ±n{x)P{x) (150) 

where Q(x) is an S0(3)-valued field, followed by over-relaxation from Q to the nearest local 
maximum qa of the DMC condition (/^ = [gAab then satisfies eq. (149)). Finally transform 
U^{x) ^U^{x) (determined by qa up to an irrelevant Z2 gauge transformation), and center 
project 

Z^{x) = signTrpt/^(x)] (151) 

to locate the P-vortices. 

All of the center gauges described here yield qualitatively similar results, but the DLC 
gauge and the method of ref. [80], which combine versions of Laplacian gauge fixing with 
relaxation to the DMC gauge, have the best center dominance properties, to be discussed 
below. 
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6 Numerical Investigations of Center Vortices 



Monte Carlo investigations of the center vortex mechanism have been mainly carried out 
in SU(2) lattice gauge theory, and the results displayed below were obtained for that gauge 
group. ^ Center projection, unless otherwise noted, is carried out in the direct Laplacian 
center gauge. 



6.1 Center Dominance 

The first question is whether center projection in DLC gauge provides a candidate for the 
decomposition 

U^{x)=U^{x)V^{x) (152) 

where U^{x) = Z^{x) (determined by center projection) carries all of the confining field 
fluctuations. The requirement is that the string tension extracted from the U = Z projected 
configuration agrees with the usual asymptotic string tension, while the asymptotic string 
tension extracted from the V configuration vanishes. We can then identify IJ (C) = ^(C*) 
as the non-confining part of loop holonomies. 

Let Z(l ^ J) denote a Wilson loop Z{C) in the projected configuration, along an / x J 
rectangular contour. The center-projected Creutz ratio is defined as 

The limit of this quantity, for large loops, is the asymptotic string tension of the center- 
projected Wilson loops Z(C). 

The results for center-projected Creutz ratios Xcp{R,R), in SU(2) lattice gauge theory 
are shown in Figs. 15 and 16. What is rather striking about these results is not only the fairly 
good agreement with the asymptotic string tension derived from unprojected loops, which 
is known as Center Dominance, but also the fact that the projected Creutz ratios are 
almost independent of R, which is known as Precocious Lineeirity. Precocious linearity 
implies that the projected potential is linear already at two lattice spacings. Another way 
of displaying both center dominance and precocious linearity is to plot the ratio 

XcpjR- R) ^ Xphys{R- R) 

as a function of the distance Rphys = Ra{(5). The result, for several values of P and lattice 
size, is shown in Fig. 17 {Xphys is Xcp in physical units, ul is in lattice units, and ^Jophys — 440 
MeV). In general, center dominance is good to about 10%. 

The infiuence of a vortex on a large Wilson loop is to contribute a factor of a center 
element, but this infiuence is only fully realized for loops whose dimensions exceed the 
thickness of a thick vortex surface. The large loops arc sensitive to the asymptotic field of 



•^Some SU(3) results in maximal center gauge are found in ref. [82], and in Laplacian center gauge in ref. 
[77]. 
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Figure 15: Creutz ratios of center-projected Wilson loops. Prom Faber et al., ref. [81]. 



vortex, which is a singular gauge transformation. Since the asymptotic effect is only seen in 
large loops, the linear potential, sensitive only to N-ality, is also only seen for large loops. 
Center projection, however, essentially collapses the thick vortex surface (or "core") onto a 
thin surface only one lattice spacing thick. This means that for the projected configuration, 
the full asymptotic effect of center vortices is felt on any distance scale greater than one 
lattice spacing. If the locations of P-plaquettes on a plane are completely uncorrelated, 
then there must be a linear potential at short distances on the projected lattice. This is the 
origin of precocious linearity. If precocious linearity is not observed, it means that either 
the P-vortex surface is very rough, or else some large fraction of the P-plaquettcs belong 
to P-vortices of small extent. Either case leads to correlations among the P-plaquettes, 
and corresponds to some high-frequency effects not directly correlated with the long-range 
physics. 

Center dominance indicates that U is a. confining background, and at least contains the 
relevant confining fiuctuations. Precocious linearity is evidence that 14 contains only the 
relevant fiuctuations. 

In the decomposition U = UV, we have assumed that the Vf^{x) fiuctuations are, by 
themselves, non-confining. This is simple to check, as first done by de Forcrand and D'Elia in 
ref. [84]. We observe that for SU(2) lattice gauge theory, with the projected configuration 
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Direct Laplacian Center Gauge 



p=2.2 
P=2.3 

W 

P=2.4 

~ m % f 

^ ^ # f I 

12* - H : 
16* ■ 

20* ' — e — I 
24* 

28* ' ' 

Michael, Teper; Bali et al. 

2 4 6 8 10 12 

R 

Figure 16: Combined data, at /5 = 2.2 — 2.5, for center-projected Creutz ratios obtained 
after direct Laplacian center gauge fixing. Horizontal bands indicate the accepted values of 
asymptotic string tensions on the unprojected lattice, with the corresponding errorbars [83]. 
Prom Faber et al., ref. [81]. 



in an adjoint gauge, 

V^{x) = Z^{x)U^{x) (155) 

V is the "vortex-removed" configuration. Figure 18 displays the Creutz ratios obtained from 
the V configuration, as compared to the center projected Creutz ratios, unprojected Creutz 
ratios, and the asymptotic string tension. Creutz ratios in the vortex-removed configurations 
tend to zero, indicating that vortex removal also removes the confinement property. 

6.2 Vortex-Limited Wilson Loops 

We define a vortcx-hmited Wilson loop Wn{C), in SU(2) lattice gauge theory, to be a 
planar Wilson loop in the fundamental representation, evaluated in the subensemble of 
configurations in which the minimal area of the loop is pierced by precisely n P-vortices (i.e. 
there are n P-plaquettes in the minimal area). Here the center projection is used only to 
select the data set; the loop itself is evaluated using unprojected link variables. 

In practice, in a Monte Carlo simulation, the procedure for evaluating Wn{C) for a rect- 
angular I X J loop is as follows: For each independent thermalized lattice configuration, fix 
to an adjoint gauge such as DMC or DLC gauge, and compute the projected configuration. 
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DLCG, Center Dominance and Precocious Linearity 
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Figure 17: Ratio of center projected and physical string tensions, vs. quark separation in 
physical units. From Faber et al., ref. [81]. 
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Figure 18: Creutz ratios in the unprojected configuration (open circles), in the center pro- 
jected configurations in DLC gauge (full circles), and in the unprojected configurations with 
vortices removed (open squares). The horizontal band indicates the accepted asymptotic 
string tension of the unprojected theory. From Faber et al., ref. [81]. 
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For each I x J loop on the projected lattice, count the number of P-plaquettes. If this 
number equals n, use the unprojected link variables along the loop to compute Tr[t/(C)], 
and include this value in the data set used to evaluate Wn{C). 

Vortex-limited loops are useful, because they allow us to check whether P- vortices on the 
projected lattice really locate the middle of thick center vortices in the unprojected lattice. 
If that is the case, then we expect in the limit of large loop area 



-1)" 



(156) 



The argument is that if center projection actually does locate correct number of vortices 
piercing a Wilson loop, then the vector potential in the neighborhood of loop C can be 
decomposed as 

A(;\x) = gdA(;\x)g-' - igd^g'' (157) 

where is a singular gauge transformation associated with the n vortices, and SA^'"^ is a 
perturbation around that background. It is understood that the delta function associated 
with the discontinuity of gdg~^ is dropped on the rhs of eq. (157). Then 



^„(C) = (-l)"(Trexp 



(158) 



It is assumed that the fluctuations 6A^^^ have only short-range correlations, and are therefore 
oblivious to the presence or absence of vortices in the middle of minimal loop area, far from 
the perimeter. In that case 



TV exp 



i ^dx^5A^^^ 



Tr exp 



i j dxi'dAf 



(159) 



for large loops, and eq. (156) follows immediately. 

Figure 19 shows our results for Wi/Wq and W2/W0, for L x L and L x {L + 1) loops 
as a function of loop area at (3 — 2.3. The data appears to agree quite well with eq. (156), 
supporting the assumption that P-vortices locate thick center vortices in the unprojected 
lattice. We also find that 



1. Xo{R, R) 0, where only the Wo{C) loops are used in computing Creutz ratios; 

2. Wodd{C) /Weven{C) ^ " Ij with loops hmitcd to having an odd or even number of 
P-vortices, respectively, piercing the loop; 

3. Xeven{R, R) 0. 



6.3 Scaling of the Vortex Density 

The asymptotic scaling property of vortex surfaces was first observed by Langfeld et al. in ref . 
[85] . Let N^or denote the total number of P- vortex plaquettes in a given lattice configuration. 
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Figure 19: Wn/Wo ratios at /3 = 2.3. From Faber et al., ref. [81]. 



and Nt the total number of plaquettes. Then the fraction of P- vortex plaquettes p is related 
to the vortex density p in physical units via 



P 



Total Vortex Area 
6 X Total Lattice Volume 



(160) 



If p is a physical quantity (i.e. constant at large /3) then the prediction of asymptotic freedom 
(inserting a(/3) into (160)) is that 

102/121 



P 



6A2 V 11 



exp 



(161) 



The average value of p can be extracted from the expectation value of center projected 
plaquettes, since it is easy to see that 

{Z{1, 1)) = (1 - p) +p X (-1) = 1 - 2p (162) 

Figure 20 is a plot of p vs. /3; the sohd line is the asymptotic freedom form (161) 
with ^ p/ 6A^ = 50. The apparent scaling of p is consistent with the density of P- vortices 
corresponding to the density of physical, surface-like objects (see also the recent results of 
Gubarev et al. [86] in the IMC gauge). Note that the scaling line has the slope appropriate 
to a density of surfaces. The scaling lines for pointlike objects (e.g. instantons) or line-like 
objects (monopoles) would have very different slopes. 
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Figure 20: Evidence of asymptotic scaling of the P-vortex surface density. The data points 
are the number of P-plaquettes per unit volume, and the solid line is the asymptotic freedom 



6.4 Finite Temperature 

The argument that center vortices lead to an area law falloff for planar loops assumes that 
the vortex piercings of the minimal area of a loop are statistically independent. It is clear 
that this assumption cannot be true for vortex surfaces whose extension is smaller than the 
size of the given loop. Let us define the extension L^ax of the vortex surface to be the 
maximum distance between any two points on the vortex sheet. In that case, a piercing of 
a plane at point a; by a vortex would be accompanied by another piercing at some distance 
L < Lmax from x. This is due to the fact that vortex surfaces are closed. If there were an 
upper limit to vortex extension (apart from that imposed by the finite size of the lattice), 
then it is not hard to show that vortex effects could at most lead to a perimeter law falloff 
for Wilson loops [87]. This means that the vortices which lead to a finite string tension must 
have an extension which is comparable to the maximum possible separation on a finite size 
lattice; we therefore expect that the vortices responsible for confinement percolate through 
the lattice. 

This brings us to the question of deconfinement at finite temperature. On time asym- 
metric lattices, the length of the lattice in the time direction plays the role of inverse temper- 
ature, and it is known that a transition from the confining to the deconfining phase in SU(2) 
lattice gauge theory occurs at a temperature around T = 220 MeV. Beyond this tempera- 
ture, Polyakov lines have a non- vanishing expectation value, and Polyakov line correlators 



prediction with Jp/QA"^ = 50. From Faber et al., ref. [81]. 
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therefore tend to a finite limit with spatial separation. If the vortex surfaces identified by 
center projection are responsible for confinement, it would make sense that these surfaces 
cease to percolate beyond the deconfinement transition. Figure 21 is a sketch of these ex- 
pectations on the projected lattice at fixed below (left hand figure) and above (right hand 
figure) the transition. In this space-slice the vortices appear as closed loops, or as loops 
closed by lattice periodicity. 





confined phase 



deconfined phase 



Figure 21: Schematic picture of P- vortices in a "space-slice" of the lattice, below and above 
the deconfinement phase transition. From Engelhardt et al., ref. [87]. 



On the other hand, even beyond the deconfinement transition, spacelike Wilson loops 
are known to retain an area-law falloff. Somehow the vortex surfaces which would disorder 

these loops must percolate through the lattice at any temperature. 

A study of P-vortex extension across the deconfinement phase transition was carried out 
by Engelhardt et al. [87] (with earlier related work carried out by Langfeld et al. [88] and 
Chernodub et al. [89]). Figure 22 displays the main result. Engelhardt et al. consider the 
extension of vortex loops in a slice of the lattice with one spatial coordinate held fixed. In 
this "space-slice" , the vortices are closed loops, and the maximum possible distance between 
links on a P-vortex loop is ^2(L/2)2 + (Lf/2)^, where L, Lf is the length of the lattice in 
the space and time directions respectively. The axis labeled "cluster extension" is in units 
of this largest distance. Figure 22 is a histogram showing the fraction of P-piaquettes in the 
space-slice belonging to vortices of a given extension. The left-hand figure is in the confined 
phase at T = 0.7 Tc, and we see that almost all vortex plaquettes belong to surfaces of 
maximal extension. The right-hand figure is at T = 1.85 Tc in the deconfined phase, and in 
this case the vortex plaquettes belong mainly to small loops. 

Spacelike Wilson loops, however, have an area law at all T. Thus, if we take a slice of 
the lattice at fixed time, we expect to see that vortices still percolate, so as to disorder these 
loops. In fact, this is just what happens, as seen in Fig. 23. Further support of this general 
picture, also obtained via center projection in DMC gauge, is found in ref. [90]. 

The deconfinement phase transition has also been investigated from the vortex free 
energy point of view by de Forcrand and von Smekal [19], Del Debbio, Di Giacomo, and 
Lucini [20], and by Hart et al. [21]. For vortex fiux in a spacelike plane (timelike vortices). 
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Figure 22: Histograms of vortex extension in a space-slice at finite temperature, both below 
(left figure, T = 0.7Tc) and above (right figure, T = f.85Tc) the deconfinement phase 
transition. The data is for = 2.4 on a 12^ x 7 (confined) and 12^ x 3 (deconfined) lattices; 
center projection in DMC gauge. From Engelhardt et al., ref. [87] 
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Figure 23: Histograms of vortex extension in a time slice at finite temperature. Vortices are 
identified by center projection at = 2.4 in DMC gauge. Data at three temperatures are 
shown on the same figure. From Engelhardt et al. ref. [87]. 



the vortex free energy goes to zero both below and above the deconfinement transition. 
However, for vortex flux in a timelike plane (spacelike vortices), the free energy goes to zero 
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Figure 24: The effect of vortex removal on tlie cliiral condensate in quenclied lattice gauge 
theory. Data points denoted "Modified" are from configurations with vortices removed in 
DMC gauge (open circles come from a looser gauge-fixing criterion). From de Forcrand and 
D'Elia, ref. [84]. 



below and infinity above the transition temperature. The interpretation of these results [19] 
is that raising the temperature by reducing the time extension of the lattice has the effect 
of "squeezing" spacelike vortices, whose flux cannot spread out in all directions of a timelike 
plane. This is what drives the deconfinement phase transition. On the other hand, timelike 
vortices are not squeezed by a reduced time extension, since their flux can spread out in a 
spacelike plane. These vortices remain condensed across the transition, and account for the 
spacelike string tension above the deconflnement temperature. 



6.5 Topological Charge 

A remarkable result, found numerically by de Forcrand and D'Elia [84], is that removal of 
vortices from lattice conflgurations removes not only the conflnement property, but chiral 
symmetry breaking disappears as well, and the no- vortex held conflgurations V^(a;) are 
always in the trivial topological sector. In Fig. 24 we show the results of these authors for 
the chiral condensate ^V-'V') a function of quark mass, at /? = 2.4 in SU(2) lattice gauge 
theory. The solid line labeled "original" corresponds to the condensate extracted from the 
Uij,{x) conflgurations, the solid fine labeled "modifled" is extracted from the corresponding 
no- vortex conflgurations V^(x) = Z*{x)Un{x) (the dashed line corresponds to a weaker 
gauge-flxing criterion). 

The fact that the topological charge of the no- vortex conflgurations V^{x) is zero suggests 
that this charge is entirely carried by the (thin) vortex background Z{C). The way in which 
topological charge can arise at isolated points on a thin vortex surface has been analyzed 
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(a) 



(b) 



Figure 25: Intersection points (a) and writhing points (b) which contribute to the topological 
charge of a P- vortex surface. Prom Reinhardt, ref. [93]. 



by Engelhardt and Reinhardt in refs. [91, 93, 75]. In general, topological charge arises at 
lattice sites at which the set of tangent vectors to the vortex surface span all four space-time 
directions. These singular points are of two varieties: 

• Self-intersections, where two areas of the vortex surface intersect at a point (Pig. 25a). 
These points carry topological charge ±|. 

• Writhing points, in which the vortex surface twists about the point in such a way as 
to produce four linearly independent tangent vectors (Pig. 25b). These points carry 
topological charge of modulus less than |. 

Supposing the center flux on the vortex surface to lie in a Cartan subalgebra of the gauge 
group, the orientation of that flux is still relevant. Vortex surfaces with net topological 
charge different from zero must be non-oricntablc, consisting of surface patches of differing 
orientation [91, 92]. The boundaries of these patches can be regarded as monopolc worldlines 
on the vortex sheet (to be discussed further in section 7, below). The net topological charge 
of a closed surface is always an integer. A very similar picture for the vortex origin of 
topological charge was put forward independently by Cornwall in ref. [94].^ 

The chiral condensate is related to the density of zero modes of the Dirac operator by the 
famous Banks-Casher formula [96] , and it is well known that quark zero modes are localized 
around instantons in an instanton background. Reinhardt et al. have shown that something 
analogous works for topological charge generated by vortices [92, 93]. The concentration of 
quark zero modes around the intersection points of four intersecting vortex sheets is shown 
in Fig. 26. 

Having associated topological charge with singular points on thin vortices, a natural 
question is whether P-vortices generated by center projection provide a topological charge 

"'In this work the boundary of a patch is identified as the worldhne of a "nexus" , which appears (along 
with center vortices) as a solution of an effective gauge theory containing a gluon mass term [95] . 
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Figure 26: Quark zero modes in a background of intersecting P-vortex sheets. From Rein- 
hardt et al., ref. [92]. 



density consistent with measurements of this quantity in the unprojected theory. Bertie, 
Engelhardt, and Faber have investigated this question in ref. [97], and find that the P- 
vortex result for topological susceptibility is consistent with conventional measurements of 
this quantity. 

Engelhardt, Reinhardt, and Faber [98] have also investigated a simple random surface 
model of vortex behavior, with a curvature-dominated action on a hypercubic lattice, which 
was put forward in ref. [99]. They report that this simple model is capable of reproducing, 
both qualitatively and quantitatively, many of the non-perturbative features of the Yang- 
Mills vacuum. Their model includes a dimensionless constant and a lattice scale, which 
are fixed by fitting to the physical string tension and deconfinement temperature. Then 
the spatial string tension at high temperature, the topological susceptibility, and the chiral 
condensate are predictions of the random surface model, and these are found to be in 
agreement with the values obtained in Monte Carlo simulations of SU(2) lattice gauge 
theory. Given the simplicity of the model, these results are of course very encouraging. A 
detailed discussion of the chiral condensate in the random surface model is presented by 
Engelhardt in ref. [100]. 

How all of this relates to the more standard scenario [101] of chiral symmetry breaking 
in the instanton liquid picture remains to be seen. Some recent data on the distribution of 
topological charge on the lattice, which may be relevant to this issue, is found in ref. [102] 
(see also [103]). 
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6.6 Casimir Scaling and Center Vortices 



The vortex mechanism is motivated by the fact that the asymptotic string tension depends 
only on the N-ahty of the color charge representation. On the other hand, as emphasized 
in section 4.2, there is an intermediate distance range, the Casimir scaling regime, in which 
string tension is instead proportional to the quadratic Casimir of the color charge represen- 
tation. 

The numerical data, however, also indicates that center vortices in SU(2) lattice gauge 
theory are comparatively thick objects, with a surface thickness on the order of one fermi. 
This estimate of thickness is based on three pieces of evidence: (i) the loop size at which 
ratios of vortex-limited Wilson loops Wi/VFq are close to —1; (ii) the spatial lattice extension 
at which the vortex free energy drops nearly to zero [16]; (iii) the adjoint string-breaking 
distance [46]. In analyzing the effect of center vortices on Wilson loops of extensions less 
than one fermi, the finite thickness of vortices must surely be taken into account. 
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Figure 27: Casimir scaling for the cooled lattice. All-loop adjoint Creutz ratios x'"^-' (i?, i?) 
(squares) compared to the corresponding Casimir-rescaled (x8/3) data for the all- loop fun- 
damental Creutz ratios (triangles). Prom Faber et al., ref. [107]. 

A very simple model of the effect of finite vortex thickness was put forward in ref. [104]. 
If the vortex core overlaps a Wilson loop perimeter, the suggestion is to represent its effect 
on the loop by insertion of a group element G somewhere in the product of link variables 
around the loop, interpolating smoothly between —I and /. Explicitly, 

G = Se-x.^[iac{x)a^/2\S'^ < aR{x) < I-k (163) 

where S is some SU(2) group element, x is the position of middle of the vortex in the plane 
of a (planar) loop C, and O-cix) — * (27r) for x in the exterior (interior) region of a large 
loop, far from the perimeter. In this simple model, S is assumed to be random on the SU(2) 
group manifold. We define / to be the probability for the middle of a vortex to pierce 
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Figure 28: Zero- vortex adjoint Creutz ratios Xq''-' (i?, i?) (squares) compared to the corre- 
sponding Casimir-rescaled (x8/3) data for the zero- vortex fundamental Creutz ratios (tri- 
angles). From Faber et al., ref. [107]. 



a given plaquette. If we then consider rectangular R x T loops (T » R) in SU(2) color 
representation j, the model predicts 



Vj{R) = - E M(i-/) + /^.Mxn)]} 

n=— oo 
1 ^ 



Gilcx] — cos(Q;m) 

■J ' m=—j 



(164) 



where Xn = n + ^ is the coordinate in the R direction. For large loops, the string tension 
has the correct asymptotic form, i.e 



{— ln(l — 2/) j — half- integer 
j — integer 

For small loops, where aR{x) << 27r, we find instead 

^.(^) = {{ E ««K)|j(j + i) 



(165) 



(166) 



which is proportional to the SU(2) quadratic Casimir (a result which generalizes to SU(N)). 
However, the potential is not necessarily linear. 

It has been found in rcfs. [104, 105] that one can choose some reasonable form of «/j(a;), 
interpolating between and 27r, such that there is an interval of charge separation R in 
which the potential is roughly linear, and roughly Casimir scaling. Without fine-tuning 
of aR{x), one can really expect no more than rough agreement from such an obviously 



60 



oversimplified model. Nevertheless, the point which is made is that the vortex theory is by 
no means incompatible with Casimir scaling, which is an efi^ect that may be attributed to 
finite vortex thickness.^ 

As a test of the relation between Casimir scaling and center vortices, one may examine 
the string tension of vortex-limited Wilson loops Wq'^-'{C) in the adjoint representation. Ex- 
cluding vortices from the middle of the loop should greatly reduce the adjoint string tension 
in the Casimir-scaling regime, if indeed vortex cores are correlated with the Casimir scaling 
effect. This test was carried out in ref. [107], which compared Xo^\R-> R) to X"''^^{R-> R)-> using 
the constrained cooling procedure [108] for noise reduction. Figure 27 shows the agreement 
with Casimir scaling for i?) in the cooled configuration. Figure 28 displays corre- 

sponding results for the zero vortex loops. Again there is Casimir scaling, but this time 
there is no apparent string tension, for either the fundamental or the adjoint representation 
loops. We conclude from this data that thick center vortices seem to be associated with 
string tension also in the Casimir scaling regime. 

6.7 Critique 

One criticism that can be levelled against the center projection data presented here is that 
the agreement of the center projected and full string tensions varies widely among different 
adjoint gauges, and among Gribov copies in the same adjoint gauge [109, 110, 111, 112, 113, 
114, 81]. DMC gauge, for example, when fixed by the over-relaxation procedure, has very 
good center dominance properties, with disagreement between full and projected string 
tensions on the order of a few percent [112]. However, using an improved procedure of 
gauge fixing based on simulated annealing, the disagreement increases to approximately 
30% [113]. An alternative to maximal center gauges is the Laplacian center gauge, which 
is free of Gribov copies. In that case the asymptotic string tensions of projected and full 
configurations seem to agree reasonably well. However, in Laplacian gauge there is no 
sign of precocious linearity or asymptotic scaling of the vortex density, and this means 
that the projected configurations, in this gauge, contain traces of irrelevant short-distance 
fluctuations. The two best methods to date, the DLC gauge and the method of ref. [80], 
begin with two different versions of the Laplacian center gauge, and subsequently "cool" 
away the irrelevant high-frequency fluctuations by relaxation to the nearest local maximum 
of Tr^[t/]. The projected configurations derived from each of these two methods, from a 
given lattice configuration, are found to be highly correlated [81]. 

Although there are some arguments about why one adjoint gauge (or set of gauge copies), 
might be better than another [77, 114, 81], the real argument that center projection locates 
vortex configurations is at present empirical; namely, the correlation of P-vortices with 
gauge-invariant observables, as seen in the results for vortex-limited Wilson loops Wn{C), 
and the observables of vortex-removed configurations — Z*!/^. 

A second area which is not entirely satisfactory is Casimir scaling. Vortex confinement 
is by no means incompatible with Casimir scaling, as discussed in the previous subsec- 

'""Wc note that the idea that the finite thickness of vortices could produce a hnear potential for adjoint- 
representation particles, over some finite distance range, was actually put forward long ago by Cornwall in 
ref. [106]. 
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tion. That is different, however, from saying that the center vortex mechanism actually 
predicts Casimir scaling, to the level of accuracy seen in numerical simulations [42, 43]. 
At the moment only crude approximations in the vortex picture, leading to equally crude 
approximations to Casimir scaling, are available. 

A third problem area concerns the Liischer 1/r term, which has not been seen in the 
center-projected data. In view of the accuracy needed to see the Liischer term in Monte Carlo 
simulations, this fact is unsurprising. Still, there is as yet no good theoretical understanding 
of the Liischer term in the context of vortex fluctuations. The effect is presumably buried in 
small correlations among vortex piercings in neighboring areas, leading to a small correction 
to the dominant area-law fall off. There does seem to be numerical evidence [115] in center- 
projected configurations, of the logarithmic growth in the width of the QCD flux tube (i.e. 
roughening), which is a definite indication of the presence of string- like fiuctuations. 

7 Monopole Confinement and the Abelian Projection 

Confinement as an effect due to abelian monopoles is one of the oldest proposals for quark 
confinement, and it is also the proposal which has received the most attention over the 
years. The idea is motivated by the squeezing of magnetic fields into flux tubes in type 
II superconductors, and by the demonstrable confinement of heavy electric charge in a 
monopole plasma, which arises in compact QED in D=3 dimensions. 

I will discuss "dual" superconductivity and monopole plasmas in a little more detail 
below. In order to apply the physics of these essentially abelian examples directly to a non- 
abelian gauge theory, it is crucial to single out an abelian subgroup, U{\)^^^, of the full 
SU (N) gauge group. This is usually accomplished by fixing the gauge so as to diagonalize 
a certain operator, leaving a U{1)^~^ subgroup. The operator may be an adjoint represen- 
tation scalar field in the Lagrangian, as in the Georgi-Glashow and Seiberg-Witten models, 
or else it may be a composite field, as suggested by 't Hooft [116] for SU(N) pure gauge 
theories. In either case, it is the electric quantum numbers of the quark with respect to 
the abelian [/(1)^~^ subgroup which are relevant for confinement in the abelian monopole 
models. 

Despite its popularity (and relative antiquity), I believe that confinement by abelian 
monopoles is problematic for QCD because it leads to the wrong represent at ion- dependence 
of the string tension. In monopole gas and dual-superconductor models, quark string tension 
is proportional to f/(l)^^^ electric charge. In contrast: 

1. At asymptotic distances, string tension depends only on N-ality. This means that for 
quarks in Yang-Mills theory with, e.g., two units of U(l) abelian electric charge , the 
string tension is zero, rather than twice the single-charge tension. 

2. At intermediate distances, the monopole prediction of U{1)^~^ electric charge depen- 
dence contradicts Casimir scaling. It also leads to an incorrect multiplicity of Regge 
trajectories. 

These remarks will be presented in more detail below, after a brief discussion of dual super- 
conductivity and compact QED^. 
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7.1 The Abelian and Dual Abelian Higgs Models 

When a type II superconductor is placed in an external magnetic field, the magnetic field can 
only penetrate the superconductor in cylindrical regions known as Abrikosov vortices, which 
carry only a certain fixed amount of magnetic flux. Abrikosov vortices are magnetic flux 
tubes, and in principle such a flux tube could begin at a monopole, of appropriate magnetic 
charge, and end on an antimonopole. The constant energy density along an Abrikosov 
vortex then implies a linear potential between the monopole and antimonopole. In other 
words, a type II superconductor, which is a system in which bosons (the Cooper pairs) of 
electric charge 2e are condensed, is a system which squeezes magnetic fields into tubes of 
quantized flux, and conflnes magnetic charge via a linear potential. 

The example of type II superconductors becomes even more suggestive in view of the 
electric/magnetic duality symmetry of Maxwell's equations. If we are willing to allow for 
the existence of elementary magnetic charges g, then the corresponding Maxwell's equations 

d^F,,=jt , a'^*F^. = jr (167) 

are symmetric with respect to the interchange of flelds and currents 

E^B , B^-E , j-^j™ (168) 

Conservation of both the electric and magnetic currents implies separate electric and mag- 
netic U(l) symmetries. It occurred independently to 't Hooft [117] and Mandelstam [118], 
in 1976, that quark conflnement in QCD might be due to some type of dual supercon- 
ductivity, in which the QCD vacuum state can be regarded as a condensate of a magneti- 
cally charged boson fleld (magnetic monopoles) , conflning electrically charged particles (the 
quarks) . 

This idea is expressed more concretely in a relativistic generalization of the Landau- 
Ginzburg theory, known as the abelian Higgs model, the action of which is 

S^Jd^x {\f^.F'^'' +\{d, + ie^^)0|'+^(00* - v'f) (169) 

The abelian Higgs model has magnetic flux-tube solutions analogous to Abrikosov vortices, 
which are known as Nielsen-Olesen vortices [119]. For a vortex along the z-axis, the 
Nielsen-Olesen vortex in cylindrical coordinates has the following asymptotic form: 

0(r)=x(r)e-^ , Ae(r) ^ A(r) 

(170) 

X(r) ''■^ V + (const.) x e ^'"^ , A{r) n/{er) + (const.) x e j yjr 



The topological stabihty of these vortices is related to the fact that the Higgs potential 
has its minimum away from = 0, and that the complex phase of the Higgs field along a 
closed loop around the vortex has a non-zero winding number. The magnetic fiux carried 
by a Nielsen-Olesen vortex is 2Tm/e. In the dual version of the abelian Higgs model, the 
Higgs fleld is magnetically charged, and it couples, not to the usual vector potential A^, but 
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rather to a "dual" photon field C^. Then the Nielsen-Olesen vortex is an electric, rather 
than magnetic flux tube, and it is electric charge which is confined. 

In both the abelian and dual abelian Higgs models, the Higgs potential has a classical 
minimum away from = 0, and the symmetry is said to be spontaneously broken. But it 
is not so obvious what that phrase really means, and which symmetry is actually broken. 
Certainly if (0) 7^ 0, then it is the local gauge symmetry which is broken, because (j) is not 
invariant under local U(l) gauge transformations. But Elitzur's theorem assures us that 
in the absence of gauge-fixing, (0) = in any phase. Of course, it is possible to fix to a 
unitary gauge so that (0) 7^ 0, but then (0) 7^ in both the broken and the Coulomb phases. 
Consider instead gauge fixing to Coulomb gauge. In that case there are some remnant global 
symmetries: a constant transformation g[x) = exp[i^], and the "displacement" symmetry 

g{x) = exp[id ■ x\ (171) 

It is breaking of these remnant global symmetries which separates the Higgs from the 
Coulomb phase in abelian theories; c.f. ref. [120] for a discussion of this point. 

There is an important limitation of the dual abelian Higgs model, which is relevant 
to its applicabihty in a non-abelian gauge context: The dual-abelian Higgs model does 
not easily accommodate fundamental fields which are electrically charged, since the kinetic 
term of an electrically- charged field involves the vector potential of the photon, not the 
dual-photon field. Despite the beautiful electric-magnetic symmetry of Maxwell's equations 
with electric and magnetic 4-currents, it is not possible to write a local Lagrangian involving 
both electrically and magnetically charged fields using only the photon, or only the dual 
photon, vector potential. And as it happens, there are electrically-charged fields in QCD 
which are highly relevant to the long-distance dynamics; these are the components of the 
gluon field which are electrically charged with respect to the U{1)^~^ subgroup. 



7.2 Confinement in Compact QED^ 

Compact electrodynamics on the lattice contains lattice-scale magnetic monopoles, with 
masses on the order of 1/e^ times the inverse lattice spacing. These massive objects interact 
with each other via a Coulombic potential. In D=4 dimensions the monopole worldlines 
form small closed loops at weak couplings, with little or no effect on the long distance 
physics, which is described by free massless photons. In D=3 dimensions the situation is 
quite different. Monopoles are instantons in that case, and the result, as shown by Polyakov 
[121], is an area law falloff for large Wilson loops. 

In the Villain version of compact QED, a series of field transformations results in the 
monopole contribution to the partition function [122] 

°° r 2n^ 1 
Zmon= Yl 5-Vm(r')G'(r-r')m(r) , (172) 

m(r)=— 00 ^ r,r' 

where m(r) is an integer valued monopole field at the (dual) lattice site r, and G{r — r') 
is the lattice Coulomb propagator (V^G(r — r') = —5rr') in three dimensions. This is 
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the partition function of a monopole Coulomb gas. The propagator can be replaced by a 
Gaussian functional integral: 



Zman = / R W exp [- ^ ^ ^ ( A^XW)^] 



°° r 27r^ 

m(r-)=-oo 9 CI r r 



(173) 



and in the weak coupling limit we need only to maintain the first terms |m(r)| < 1 in the 
sum, so that 



7 ~ 
^mon ~ 



\{dx{r) exp[-|-^^(^(A^x)^ - cosx(? 



47r2 V^2' 



where 



Mo' = — exp 



G(0) 



(174) 



(175) 



Going to continuum notation, the action is 



/ ^x(r) 



exp 



47r2 



/ d\ {^{d^xf - M2cosx(r))]. (176) 



where M = Mo/a. 

The electromagnetic field due to the monopole density is given by 



d^H^ir) = 27rm(r), i.e. H^{r) ^^J dh' 



./|3 



m r 



:i77) 



Then if C denotes a closed curve and S{C) a surface with C as boundary, we have that 

/■ dr^A^ir) = / dS^ir) H^{r) = / dh r^sic){r) m(r), (178) 



where 



^5(c)(r) 



1 a 



/ rf5^(r') ^— . 
JSiC) \r — r'\ 



(179) 



2 i5(C) 

The monopole contribution to a Wilson loop carrying n units of the elementary electric 
charge is 

Un{C) = e*"/*> = / ^te(c)('') ™W (180) 
Inserting this expression into the partition function for the monopole Coulomb gas (173) 

-.2 



mc)) 



Vx{r) exp[-|^ I d^r (^^{df,{x - nr]s(c)f - cos x{i 



[181) 



The functional integral is then evaluated by a saddlepoint method, solving the equations 
of motion for x in the presence of a source nr]s(c)- The calculation for n = 1 was first carried 
out by Polyakov [121, 122] and the result is 



(C/i(C)) exp —a area(C) 
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(182) 



where area(C) is the minimal area of the loop. The Polyakov result has a simple general- 
ization to arbitrary integer charges, found in ref. [123], 

(183) 

which will be of some importance in our discussion of the monopole mechanism in non- 
abelian gauge theories. We note that compact QED can also be regarded as a particular 
limit of the dual abehan Higgs model [124, 125]. 

For compact QED in D = 4 dimensions, the integer-valued monopole field m(r) in 
the monopole action (172) is replaced by an integer- valued, conserved monopole current 
A;^(r), defined on links of the dual lattice. Smit and van der Sijs [125] proposed that an 
abelian monopole action of this type, which in addition includes a local monopole mass term, 
could be regarded as the effective long-range action of a pure non-abelian gauge theory 
in D=4 dimensions. This effective action is deduced from the action of multi-monopole 
configurations in an underlying Yang-Mills theory. In the absence of a Higgs field, BPS- 
like monopole configurations [126] in the pure Yang-Mills theory were considered. In this 
approach the concept of a monopole "container," i.e. a finite volume with certain monopole 
boundary conditions, is introduced to specify the relevant configurations in the quantized 
theory, and to compute some of their properties [125, 127]. 



{Un{C)) exp —na area(C) 



7.3 The D=3 Georgi-Glashow Model 

The Georgi-Glashow model, aka the SU(2) adjoint Higgs model, is an SU(2) gauge theory 
with a Higgs field in the adjoint representation. In D = 3 dimensions the action is 



S = jd'x [-Tr [Fl^] + -{D,4>''f + -^Krr - v'f 



(184) 



The model has charged W-bosons of mass Mw = gv, a massive Higgs field of mass Mh = 
-\/2Av^, a (perturbatively) massless "photon," associated with the unbroken U(l) symmetry, 
and 't Hooft- Polyakov monopoles (which are instantons in D=3 dimensions), with action 



^2 VMi 



LW 



where e(r) is a slowly varying function of order one. In a unitary gauge, where the Higgs 
field is rotated to point in the positive color 3 direction, the field is the photon, and the 
charged W-bosons are the — A^^± lA^^ fields. In unitary gauge, the 't Hooft-Polyakov 
monopole has the form of a Dirac monopole in the A^^ field component 

Fl^e,^^d^Al^-^ (186) 

at distances far from the monopole center (where |(/)| = v). The 't Hooft-Polyakov monopoles 
again interact via a Coulombic potential, and the corresponding partition function has the 
form of eq. (176), except that now the parameter M, which diverges in compact QED^ in 
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the lattice spacing a — limit, is replaced by a finite constant, dependent on the coupling 
g and W-mass Mw Wilson loops can be computed as in compact QED^, and a finite 
string tension a ~ exp[— ^^o^] is obtained. It would seem that in the Georgi-Glashow 
model, the goal of separating the gauge field into a piece A, which carries only the confining 
fiuctuations, and the remaining non-confining fiuctuations A = A— A has been accomplished 
unambiguously: The instruction is to go to unitary gauge, identify the monopole positions 
{x'^} from the zeros of the Higgs fields, and let «4^(a;) = S^'^A/^ be the vector potential of 
Dirac monopoles at those positions, so that 



Un{C) = exp 



A 

where Qk — il/g'- This identification of the confining fiuctuations is certainly correct for 
compact QED^. It is not correct, however, in the Georgi-Glashow model. 

The problem is connected with loops of multiple electric charge. Consider, in unitary 
gauge, loops of the form 

in- i Al{x)dx^' (188) 

If the confining fields arc those generated by a monopole Coulomb gas, then the charge- 
dependence of the string tension is the same as for compact QED^, namely 

(UniC)) ~ e-'^'^^(^) (189) 

with A{C) the minimal area. But this answer cannot be correct in the Georgi-Glashow 
model, as pointed out in ref. [123], because there are charged bosons carrying two units 
of the minimum electric charge. This means that static charges carrying n units of electric 
charge are screened to n mod 2 units of charge, and therefore 

The vector potential obtained from the monopole Coulomb gas action leads to the wrong 
answer for these observables, and from this it follows that confining magnetic flux Al^{x) is 
not distributed in accordance with a monopole plasma. 

What we learn from this is that the infrared effects of W-bosons cannot be neglected 
in constructing an effective action for the Georgi-Glashow theory. The usual justification 
for ignoring the W-bosons is that these are massive objects, and therefore cannot affect the 
dynamics at large scales. For a confining theory, that argument is simply wrong. Moreover, 
confining fiuctuations must organize themselves at large scales in such a way that n — even 
charged loops have a vanishing string tension. The most obvious way to satisfy this condition 
is for the confining flux to be coUimated into Z2 vortices, rather than being distributed as 
in a monopole Coulomb gas [123]. 

On the other hand, the monopole Coulomb gas effective action (176) ought to have 
some range of validity in the D=3 Georgi-Glashow model. On the basis of simple energetics 
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arguments, the electric flux tube between static sources carrying two units of U(l) electric 
charge will break by W pair creation only when the energy of the flux tube is greater 
than the combined mass of the two W bosons; i.e. when La2 > 2mw The double-charge 
string tension (72 is zero beyond the string-breaking distance; however, up to that point, the 
monopole Coulomb gas result that (T2 = 2ai should hold. Then the string-breaking distance 
is L ^ m]v/(7i, and up to that length scale a monopole Coulomb gas analysis ought to give 
a good account of the dynamics. 

For the D=3 Georgi-Glashow model there exists a proposal for the infrared effective 
theory, in terms of elementary vortex field operators, which reduces to the monopole action 
(176) in a certain approximation. W-bosons, in this dual formulation, appear as solitons of 
the effective action. This effective theory was originally proposed by 't Hooft in ref. [13], and 
has been investigated extensively by Kovner and co-workers [128]. The action is motivated 
by 't Hooft's observation that the analog of the SU(N) vortex creation operator -B(C), in 
2-1-1 dimensions, is an operator V{x) defined such that 

Vix)\AkA) = \Al<l)') (191) 

where A', 0' are related to A, (p via a singular gauge transformation, defined such that 

W{C) e±2WA'W^(C) (192) 

if C is topologically linked to point x (with the ± sign depending on the clockwise or 
counterclockwise orientation of the loop). In the Georgi-Glashow model, suitably generalized 
to an SU(N) gauge group, V{x) creates a topologically stable soliton, and the number 
of solitons minus antisohtons is conserved mod N (the worldline of a such a soliton is a 
center vortex in D=3 dimensions). This means that SU(N) adjoint Higgs theory in D — 3 
dimensions has a global Z^v topological symmetry associated with this conservation law, 
and 't Hooft argued that confinement is associated with the spontaneous breaking of this 
dual topological symmetry.^. On these grounds, 't Hooft proposed an effective action for 
the soliton field with a Lagrangian 

L = d^Vd'^V* - X(V*V - + m(v^ + {V*f) (193) 

Kovner ct al. [128] point out that L can be related to Polyakov's effective action in eq. (176) 
by taking N — 2 and A very large, so that the modulus of V is almost frozen. Writing 



V(x) — nexp 



(194) 



and inserting this expression into eq. (193) yields the Polyakov action in (176), provided 
one ignores the fact that x(a;) is actually a phase with a certain periodicity, rather than a 
continuous field. In the 't Hooft action (193), as opposed to the Polyakov action, the field 
x{x) is allowed to have quantized discontinuities, which are essentially vortex configurations 
of the vortex field itself. It turns out that these vortices in the vortex field arc associated 
with charged W-bosons, whose importance to infrared dynamics has already been noted. 
Further details concerning this vortex operator approach to confinement in the D = 3 
Georgi-Glashow model can be found in the review articles by Kovner et al. [128, 129]. 

^Of course the usual Zn global center symmetry, discussed in section 4.2, is unbroken in this phase 
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7.4 The Abelian Projection 



In the Georgi-Glashow model, the abehan U(l) subgroup is singled out by the adjoint 
Higgs field. In the unitary gauge 0" = p5"^, the "photon" of the theory is simply the 
component of the gauge field, and U(l) magnetic monopoles are associated with singularities 
in the gauge-fixing condition; i.e. points (D=3) or lines (D=4) where = 0. In QCD, on the 
other hand, there are no scalar fields in the Lagrangian, and no obvious way of distinguishing 
an abelian subgroup. 

In ref. [116] 't Hooft suggested that the adjoint scalar field used to identify the abelian 
subgroup could be a composite operator X, formed from the gauge field. Operators such as 
X = F^j^ or Fi2, which transform under a gauge transformation g via 

X{x) ^ X'{x) = gix)X{x)g-\x) (195) 



are candidates for this role. The corresponding unitary gauge, in an SU(N) gauge theory, 
is the gauge in which X is diagonalized, i.e. 



Ai 



A. 



X = 



(196) 



A 



N 



and the eigenvalues ordered Ai < A2 < ... < Ajv- This gauge choice leaves a remnant U{1)^~^ 
gauge symmetry, generated by the Cartan subalgebra of the SU(N) gauge group. The 
"photons" corresponding to the remnant abelian gauge symmetry are the diagonal elements 
A'Jf'^x) of the matrix- valued SU(N) vector potential, while the "W-bosons" correspond to the 
off-diagonal elements. Singularities of the gauge-fixing condition are identified with points 
(D=3) or fines (D=4) where two of the eigenvalues Aj of X coincide; these singularities are 
identified as magnetic monopoles of the abelian subgroup. All this is in complete analogy to 
the Georgi-Glashow model, except that a composite operator X replaces the adjoint Higgs 
field. As in the dual superconductor picture and in compact QED, confinement would be 
associated with the condensation of monopoles, and the confining fiuctuations that we have 
denoted by An would be entirely contained in the diagonal ( "photon" ) components of the 
vector potential. 

In translating this idea to the lattice formulation, and testing it by Monte Carlo sim- 
ulations, the first issue is the choice of gauge. In 't Hooft's original article, the emphasis 
was on "non-propagating" gauge choices, which are free of massless ghosts. Within this 
class, any gauge choice would suffice. This gauge independence does not appear to survive 
lattice tests, where strong correlations have been found between, e.g., the gauge choice and 
the string tension arising from the diagonal ( "photon" ) piece of the vector potential. Most 
lattice investigations have been carried out for the SU(2) gauge group, and we will now focus 
on this particular case. The most widely used and successful gauge choice in SU(2) lattice 
gauge theory is the McLximal Abelian gauge [76] , defined as the gauge which maximizes 
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the quantity 



D 



i? = $:ElY[C/^(x)<73C/t(x)<73] 

X fj,= l 



(197) 



The trace can be recognized as the 33 component of the U^^{x) hnk variable in the adjoint 
representation, and the condition that R is maximized is simply the condition that the link 
variables are, on average, as diagonal as possible. As with maximal center gauge, there is 
no known method for finding the global maximum of R. Instead one obtains local maxima 
(Gribov copies) from a combination of simulated annealing and over-relaxation techniques. 
Locally, in any Gribov copy of maximal abelian gauge, the operator 



y{^) = E [u^{x)a'ui{x) + ui{x - n)<j%{x - n) 



^l=\ 



(198) 



is diagonal. In the continuum, the diagonality condition is 

{d, ± i^4)Aj = (199) 
On the lattice, maximal abelian gauge allows the residual U(l) gauge transformations 

U^{x) e^'^(^)'^3[/^(^)g-i0(x+MV3 (200) 

Define u^{x) as the diagonal part of the link variable in maximal center gauge, rescaled so 
as to restore unitarity. Writing SU(2) hnk variables in terms of Pauh matrices in the usual 
way 

U^{x) ^ a^I + ia- B (201) 

this means 



Uf,{X} 



je^^ix) 



Q-i9)i{x) 



(202) 



The full link variable U^{x) can be factored into the diagonal matrix u^j, and another matrix 
Cn which contains off-diagonal elements 



U^{x) = Cfj_Uf,{x) 



(1- 



1/2 



Cfiix) 

(i-Icm(^)P 



1/2 



(x) 



(203) 



It is not hard to see that under the remnant gauge transformation, the diagonal link variables 
Ufj_ transform like an abelian gauge field 



Unix) 



^i(t>(x)(J3 



u^{x)e 



i(f){x+li) 



T3 



(204) 



70 



while c^(x) transforms like a matter field with two units of electric charge 

c^{x) ^ e2^'^(^)c^(x) (205) 

After maximal abelian gauge fixing, the mapping from the original SU(2) link variables to 
the U(l) hnk variables 

U^{x) ^ u^{x) (206) 

is known as the "Abelian Projection" . Wilson loops constructed from the abelian pro- 
jected hnk variables u^{x) are referred to as "abelian projected" (or simply abelian) Wilson 
loops. 

Abelian Dominance is the property that the asymptotic string tension extracted from 
abelian projected Wilson loops agrees with usual asymptotic string tension. This property 
was first demonstrated by Suzuki and Yotsuyanagi [130]; the most accurate measurements 
to date find an abelian projection string tension at (92 ±4)% of the full SU(2) string tension 
[131]. 

Abelian dominance is not yet proof of a monopole confinement mechanism; the confining 
fiuctuations could still, e.g., take the form of vortex configurations.^ Nor docs the abelian 
projection, by itself, separate the (infrared) confining from the (ultraviolet) perturbative 
field fiuctuations, since the abelian projected potential at short distances, which is due 
presumably to onc-photon exchange, is not linear. 

An attempt to go beyond abelian projection, and to really isolate the confining back- 
ground in the abelian-projected lattice configuration, was initiated by Shiba and Suzuki in 
ref. [133], and by Stack, Nieman, and Wensley in ref. [134]. These authors begin by iso- 
lating monopoles in the abelian-projected lattice configuration according to the De Grand- 
Toussaint criterion. Let 

fnu{x) = d^e^{x) - dj^{x) 

= 7^,(x)+27rn^,(x) (207) 

where 9^{x) is related to the abelian link variables according to eq. (202), is the forward 
lattice derivative, f ^y{x) is in the range [— tt, tt], and 71^^ is an integer valued "monopole 
string" variable. The monopole current on the lattice, according to De Grand and Toussaint 
[135], is given by 

k^{x) = —e^ap^djf^{x) (208) 
47r 

The next step is to construct the lattice configuration corresponding to a monopole 
Coulomb gas, with the monopole sources specified by the current kn{x) above. This config- 
uration is determined by the monopole string variables n^,^ as follows 

«7"(a;) = exp[^^^7"(x)] 

^r"(^) = -J2Dix-y)dln^,iy) (209) 

y 

^Actually, even a "naive" abelian projection, carried out with no prior gauge-fixing whatever, has the 
property of abelian dominance. This is for rather trivial group- theoretic reasons, however, which seem to 
have nothing to do with the confinement mechanism [132]. 
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where &^ is the backward lattice derivative and D{x) is the lattice Coulomb propagator. 
The "photon" contribution is then defined as the difference between the abelian projected 
and monopole configurations 

ul\x) = exp[z^f (x)] 

ef{x) = e,{x)-e^'^-{x) (210) 
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Figure 29: Static potentials computed from monopole fields only and photon fields 

{yreg\^ Only, in units of lattice spacing a — 0.081 fm. The potential computed from the 
abelian projected lattice is denoted V^^. Prom Bah, ref. [142]. 



Figure 29 shows the static potential derived from abelian projected configurations (V^^), 
the monopole Coulomb gas configuration (y*'"^), and the photon "remainder" field {V^^^) 
[131]. This figure certainly suggests that the monopole Coulomb gas picture of Shiba and 
Suzuki is on the right track. The monopole potential is linear starting from one lattice spac- 
ing, which means there are no extraneous high-frequency fluctuations, and the associated 
string tension is about 95% of the abehan projection string tension (itself about 92% of 
the full SU(2) string tension). This close agreement is known as Monopole Dominance. 
The Wilson loops constructed from u^^ alone seem to have no string tension at all, which 
would appear to confirm their role as non-confining, mainly high-frequency, fluctuations. 
Although total monopole density doesn't scale, the density of monopoles belonging to the 
largest cluster of monopole worldhnes does scale according to asymptotic freedom [136]. 
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Moreover, the position of monopoles identified by the abehan projection is highly correlated 
with the gauge-invariant action density of the full SU(2) gauge fields at those locations [137]. 

Further investigations have gone in several directions. Suzuki and co-workers have con- 
centrated on constructing effective actions in terms of monopole currents and their Coulom- 
bic interactions, which reproduce the monopole dominance results obtained from abelian- 
projected configurations [138]. There have also been detailed comparisons of the electric flux 
tube, found in abelian- projected lattice configurations derived from lattice Monte Carlo sim- 
ulations, with the form which is expected in the dual abelian Higgs models [139, 140, 141]. 
These comparisons work out rather well. For details, the reader is referred to the reviews 
in refs. [142, 143, 124] 



7.5 Monopole Operators 

There has been some effort, largely by the Pisa group [144], to verify dual superconductivity 
by demonstrating that a suitably defined monopole creation operator, denoted /i, acquires 
an expectation value in the confinement phase of gauge theories at finite temperatures. 
The idea is that if is magnetically charged, and acquires a VEV, then a magnetic U(l) 
symmetry is spontaneously broken. 

The n{x, t) operator is defined so as to shift any given gauge field configuration Ai^{y, t) 
at a given time t by the field bk{y;x) of an abelian monopole located at point x. The 
U{1)^~^ subgroup of SU(N), required to define the monopole, is chosen according to some 
abelian projection. In practice, what is calculated is 

P=-^H^I')]=i^)s-{S + ^S)s^^s (211) 

where S + AS is the action obtained by shifting the link variables by a monopole field at 
time t. {S)g is the usual plaquette action, and {S + A.S)g_^_^g is the expectation value of 
the shifted action, where the shifted action is also used for the Boltzmann distribution. The 
calculation has been carried out for lattice gauge theories at finite temperatures, and it is 
found that the transition from (fi) 7^ to (fi) = indeed coincides with the deconfinement 
transition, for gauge theories with no matter fields.® Related results were obtained by Cea 
and Cosmai [145]. 

Since the fi operator makes no obvious reference to the center of the gauge group, a 
natural question is what happens when matter fields are added to break the global center 
symmetry explicitly. Let us again consider the Fradkin-Shenker model. If (fi) = corre- 
sponds to deconfinement, then, since there is no Higgs/confinement transition in this model, 
one would expect that (n) — throughout the phase diagram. Unfortunately no results 
have yet been reported for the Fradkin-Shenker theory. There are new results, however, 
obtained recently by the Pisa group, for (fi) in QCD with dynamical quarks [146]. It is 
found that (/i) 7^ at low temperatures, which is contrary to what ought be the case in the 
Fradkin-Shenker model. A transition is also found to (/i) = at finite temperature in full 



^Actually, the same result is found when fj, is defined to add a vortex configuration, rather than a 
monopole field, c.f. ref. [20]. 
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QCD, in a region of the phase diagram where it is generally believed that only crossover 
behavior, rather than a true transition occurs. It is probably fair to say that the situation 
is not completely clear at the present time.^ 

7.6 Objections 

As discussed at some length in section 4, string tension varies with the group representation 
of the color charges. At asymptotic distances, the string tension depends only on the N- 
ality of the representation, while in an intermediate regime of charge separation, up to 
the onset of color screening, the string tension is proportional to the quadratic Casimir of 
the representation. These dependencies are different from what is expected (and found) in 
monopole Coulomb gas and dual abelian Higgs models, where the string tension is simply 
proportional to the abelian electric charge relative to the U{1)^~^ gauge group. Moreover, 
the monopole and dual-superconductor models predict a multiplicity of different types of 
electric flux tubes, in the intermediate distance regime, which should not exist in QCD. We 
will consider the intermediate and asymptotic distance scales separately. 

7.6.1 Intermediate Distances 

As shown in ref. [123], the string tension of a charge- anticharge pair in compact QED3, 
where the charge of each source is n units of the minimum charge, is simply 

CTn = nai (212) 

The result can be related to a picture of n independent electric flux tubes running between 
the two charges. A similar result is expected in U(l) monopole Coulomb gas models in 
D=4 dimensions, and in dual (and real) superconductors. Numerical simulations of double- 
charged abelian-projected Wilson loops, i.e. 

^ {exp[2ie{C)]) (213) 

do in fact display string tensions which are roughly double the single-charge string tension 
[149, 142], where 9{C) is the lattice "loop integral" of link angles 9^{x) around loop C. 

The problem is that if the conflning force is sensitive only to abehan electric charge, 
then it is difficult to understand Casimir scaling, and in particular how there could ever be 
a linear static potential between color charges in the adjoint representation [35]. Quarks in 
the adjoint representation of the SU(2) color group come in three different varieties, with 
+2, 0, — 2 units of the minimal U(l) electric charge, respectively. The zero electric charge 
component should completely dominate the behavior of the adjoint Wilson loop, and these, 
according to dual superconductor and monopole gas ideas, are impervious to the conflning 

^Frohlich and Marchetti have pointed out that the Pisa monopole operator actually depends on the 

position of the associated Dirac string [147]. They have suggested a related construction for the monopole 
operator which avoids this problem, and which makes use essential use of center vortices. Some numerical 
work, based on this construction, has been reported in ref. [148]. 
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force. In the abelian projection 



1 



(Tr [Ua{C)]) 



^ (exp[2im^(C)]) 



m=— 1 




,2ie(C)\ , /„-2i9(C) 



(214) 



The zero-charge component contributes a constant (one) to the trace in the adjoint repre- 
sentation. The neglected off-diagonal gluons would presumably convert the constant to a 
perimeter-law falloff, but in any case there is no reason to expect an area-law falloff for zero 
N-ality Wilson loops, much less Casimir scaling. 

A second problem at intermediate distance scales, pointed out by Douglas and Shenker 
in ref. [44], is that in the SU(N) case the corresponding abelian projection is to a U{1)^~^ 
theory, which in the dual superconductor picture would involve — 1 copies of the dual 
abelian Higgs model, each with its own solitonic flux tube. The number of flux tubes of each 
type would be a conserved quantity, which leads in principle to physically distinct Regge 
trajectories which are not expected in SU(N) gauge theory, and not observed in QCD. 

7.6.2 Asymptotic Distances 

The problem at asymptotic distances is that the string tension 02 of double-charged abehan 
Wilson loops should be zero, not 2ai. This is because, in addition to monopoles and 
"photons," the theory contains off-diagonal gluons carrying two units of the minimal U(l) 
electric charge. Eventually it becomes energetically favorable for the off-diagonal gluons 
to screen the static double charges, and the string tension beyond that screening length 
vanishes. In general, for the U(l) projection of SU(2) gauge theory, the asymptotic string 
tension between sources carrying n units of U(l) charge is 



It follows that the abelian fleld distributions characteristic of a dual superconductor, or a 
monopole Coulomb gas, cannot be the distribution of conflning flelds in the abelian projec- 
tion. 

It is difficult to check double charge screening numerically; this has only recently been 
done for SU(2) adjoint Wilson loops [48], as already noted. Instead, the expectation values 
of double-charged Polyakov lines have been computed for abelian projected configurations, 
and in the corresponding monopole dominance configurations (i.e. fields obtained from eq. 
(209), from monopole string variables identified in the abehan projection). The abehan 
projection and monopole dominance results for the double-charged Polyakov line are dras- 
tically different, as shown in Fig. 30. The data shown in this figure were obtained on a 
12^ X 4 lattice, at various values of 13. The deconfinement phase transition, on a lattice of 
four lattice spacings in the time direction, occurs around (3 = 2.2. The monopole dominance 
data is consistent (or nearly so) with a vanishing double-charged Polyakov line, implying 




n even 



n odd 



(215) 
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confinement of quarks with two units of electric charge. In the abehan-projected configu- 
rations, in which the abehan field distribution is determined by all the degrees of freedom, 
including the off-diagonal gluons, it is very clear that the double-charged Polyakov line has 
a non-zero expectation value in the confined phase; i.e. double electric charges are screened 
rather than confined. 
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Figure 30: Double-charged Polyakov lines, computed in the abelian projection and in the 
monopole dominance approximation, on a 12^ x 4 lattice. From Ambj0rn et al., ref. [150]. 



What is going wrong, in the monopole Coulomb gas and dual abehan Higgs models, is 
that they neglect the long-range effects of the off-diagonal gluons (the "W-bosons"). It is 
assumed in these models that the off-diagonal bosons acquire a mass, and therefore their 
effects on the abelian gauge field can be neglected at large scales. This is clearly incorrect. 
So what would happen if the off-diagonal gluons are properly included? A confining theory 
which contains abelian magnetic monopoles and dynamical electric charges is neither a 
monopole Coulomb gas, nor a simple dual superconductor at large scales.^" In the case 
of Yang-Mills theory in the abelian projection formulation, the vacuum distribution of the 
abelian projected field is determined by an effective action 



exp 



Seff[A^^ = j D A (ghosts) e-^^^-^«^ (216) 



so that the vacuum distribution of the "photon" field is determined, in part, by virtual 
"W-bosons" A^. This is in contrast to the monopole dominance approximation (209), where 
only the distribution of monopoles, but not the distribution of the fields emanating from 
those monopoles, is affected by the W-bosons. In general, it is clear that the W-bosons must 



^°In particular, there are no local Lagrangians which are able to couple an abelian gauge field to both 
electrically charged and magnetically charged matter fields. 



76 



have a back reaction on the vacuum distribution of the "photon" field A^^, in such a way 
that the vacuum fluctuations of the fleld are consistent with the Z2 dependence shown 
in eq. (215). 

The argument that the confining abehan configurations generated by Seff[A^] must be 
vortex configurations is the same as in section 5. We assume that vacuum configurations 
can be regarded as the sum of a high-frequency fluctuation Aj^, responsible for most of the 
perimeter-law behavior, on a conflning background .4^, which is mainly responsible for the 
area-law behavior. Then, since a2n — 0, it must be that the confining fiuctuations have 
almost no effect on the expectation values of even powers of the abelian loop holonomy 
U{C) — exp[i f which is accomplished if 



[/'"(C) ^ f/'"(C) 



C/(C) = Z(C)f/(C) 



(217) 



where in this case Z{C) = ±1. So the confining fluctuations A^^ arc those for which W(C) = 
Z{C) for every large loop. These, of course, are Z2 center vortices, and on the lattice 
correspond to link variables U,j,{x) — ± exp[i(0(x) — (/)(x -|- /i)]. 

An instructive example is the case of compact QED4 at strong coupling. The theory can 
be rewritten in terms of monopole currents interacting via a Coulombic potential (as in the 
D=3 case discussed in section 7.2), and conflncmcnt can again be understood in terms of a 
monopole Coulomb gas. However, the situation is changed quaUtatively by the addition of 
a charge- 2 matter field. For simphcity we consider a charged scalar matter field p of fixed 
modulus IpI = 1, and 



Z^ [ DpD9^ exp /35^cos(^(p)) + ^A^{p*(x)e''^''(^V(2^ + Ai) + c.c} 



(218) 



with P ^ 1 (confinement) and A ^ 1. Integrating over the matter field generates a set of 
charge-2 Wilson loops in the induced action, which screen any external even-charged Wilson 
loop. 

This gauge-|-matter theory can be rewritten, following Chernodub and Suzuki [151], in 
terms of monopole currents interacting with charge-2 electric currents induced by the matter 
field. Again, it is the electric currents which screen an external charge-2 Wilson loop. In this 
case, however, rewriting the theory in monopole variables actually obscures the underlying 
physics. The confining field configurations are no longer Coulombic fields emanating from 
monopole charges. Rather, the confining configurations are thin Z2 vortices — a fact which 
is invisible in the monopole formulation. To see this, go to the unitary gauge p — 1, which 
preserves a residual Z2 gauge invariance, and make the field decomposition 



where 
and 



exp[i9^{x)] = Zi^{x)e^[i9i^{x)] 
z^(x) = sign[cos(6'^(x))] 



n E 

^,1^ Zij,{x)=±l 



tt/2 



27r 



exp 



P Z{p) cos{9{p)) + A ^ cos(2^^(x)) 



(219) 
(220) 

(221) 
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This decomposition separates lattice configurations into vortex degrees of freedom (the 
Zp,{x)\ and small non-confining fiuctuations around these vortex configurations, strongly 
peaked at ^ = 0. One can easily show, for /5 -C 1, A 1, that 



;exp[m^(C)]) « (^"(C)) (exp[me(C)]) (222) 

with 



exp[—(T 74(C)] n odd 
1 n even 



(Z-(C)) = 

exp[m6'(C)]\ = exp[-/xn'P(C)] (223) 



which establishes that the confining fiuctuations, in this coupling range, are entirely due 
to thin vortices identified by the decomposition (219) in unitary gauge. It is clear that 
the addition of a charge-2 matter field has resulted in a qualitative change in the physics 
of confinement. Yet the transition from a monopole Coulomb gas mechanism to a vortex 
dominance mechanism is essentially invisible if the gauge+matter theory is rewritten in 
terms of monopole + electric current variables [151]. which in this case tend to obscure, 
rather than reveal, the nature of the confining fluctuations. 
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Figure 31: Monopole confinement models in SU(N) gauge theories predict that confine- 
ment forces are proportional to U{\)^~^ electric charges. Center vortex theories predict a 
dependence only on N-ality. 



For a variety of theories, the dependence of string tension on group representation is 
summarized in Fig. 31. In compact QED^ with no dynamical matter fields, the static 
potential is perturbative at short distances, followed by a confining region where the string 
tension is roughly proportional to the electric charge. This dependence of string tension on 
electric charge, which is characteristic of monopole Coulomb gas mechanisms, also applies 
to the dual abelian Higgs model. In the D=3 Georgi-Glashow model and Seiberg-Witten 
models (see below), which both have matter fields in the adjoint representation, there is 
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an intermediate distance range in which string tension is proportional to U{1)^~^ electric 
charge, followed by an asymptotic regime where string tension depends on N-ality. Finally, 
in QCD, there is no U{1)^^-^ region at all. If there were, then adjoint representation Wilson 
loops would be dominated by their C/(l)^~^ charge neutral component, and would have zero 
string tension at all distance scales. The transition is from Casimir scaling at intermediate 
distances to N-ality dependence at large distances. In QCD, there is no distance interval 
where the representation-dependence of string tensions can be directly attributed to a dual 
abelian Higgs model, or a monopole Coulomb gas. 

7.7 The Seiberg-Witten Model 

The important work of Seiberg and Witten on duality and confinement in supersymmetric 
gauge theories is beyond the scope of this review, and the reader is referred to the original 
article [152] and also to ref. [153] for details. A few comments about their work, however, 
are relevant at this point [150]. 

M = 2 super Yang-Mills theory, like the Ceorgi-Glashow model, has a scalar field in the 
adjoint representation, which can be used to single out a compact abelian subgroup: U(l) 
for super Yang-Mills theory, U{1)^~^ for the SU(N) gauge group. Seiberg and Witten were 
able to show that there is a point in moduli space at which the magnetic monopoles of this 
theory become masslcss. On adding a soft supersymmetry breaking term, which reduces the 
M = 2 supersymmetry to A/" = 1, the theory goes into a confining phase, and this transition 
is associated with a condensation of the monopole field. Seiberg and Witten derive in this 
case an effective low energy action which is a supersymmetric generalization of the dual 
abelian Higgs model. 

However, just as in Polyakov's monopole gas treatment of the D=3 dimensional Georgi- 
Glashow model, the Seiberg-Witten low-energy effective action neglects the W-bosons of 
the theory, on the grounds that they are massive and should not contribute to long-range 
physics. But we have already seen that the W-bosons are very relevant to long-distance 
physics; the N-ality dependence of string tensions Ur cannot be obtained without them. 
This means that the Seiberg-Witten effective action is not the whole story at large scales. 
This effective action is obtained keeping only local terms with no more than two derivatives 
of the bosonic fields. It is not obtained by actually integrating out the massive W-bosons, 
for otherwise the effective action would have the correct N-ality dependence built in. At 
distance scales on the order of the color screening length, L ^ mw/cr-, the dual abelian Higgs 
action arrives at the wrong N-ality dependence for Wilson loops, and is simply inadequate 
to describe the large-scale vacuum fiuctuations. This situation is essentially as depicted in 
Fig. 31: there is an intermediate range of distances, up to the screening scale, for which the 
Seiberg-Witten effective action gives a good account of the physics. Beyond that scale, the 
contributions of virtual W-bosons to the effective action can no longer be ignored. 

7.8 Monopoles and Vortices 

Let us now consider how center vortices would appear in the abelian projection (cf. refs. 
[154, 150], and also ref. [77]). We begin with a center vortex at a fixed time, as indicated 
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schematically in Fig. 32. In the absence of gauge fixing, the vortex field strength points 
in arbitrary directions in color space. Upon fixing to maximal abelian gauge, the vortex 
field strength tends to line up, in color space, in the icr^ direction. But there are still 
going to be regions along the vortex where the field rotates from the to the —a^ 
direction, as shown in Fig. 33. Upon abelian projection, these transition regions show up as 
monopoles or antimonopole (Fig. 34). The end result is that a center vortex appears, after 
abelian projection, as a monopole-antimonopole chain, with a typical vacuum configuration 
indicated very schematically in Fig. 35. If this picture is correct, then the ±27r monopole 
flux is not distributed symmetrically on the abelian-projected lattice, as one might expect 
in a Coulomb gas. Rather, it is coUimated in units of ivr along the vortex line. Because 
of this collimation, monopole magnetic fields will not affect double-charged abelian Wilson 
loop (or any abelian loop with an even number of units of electric charge). 

Figure 32: Vortex field strength before gauge fixing. The arrows indicate direction in color 
space. 



Figure 33: Vortex field strength after maximal abelian gauge fixing. Vortex strength is 
mainly in the ±(73 direction. 




Figure 34: Vortex field after abelian projection. 

It is important to see whether this collimation of field strength along the vortex line, and 
the tendency of abelian monopoles to line up along vortex lines in a monopole-antimonopole 
chain, can actually be observed in Monte Carlo simulations. This question was studied in 
refs. [154, 150] in SU(2) lattice gauge theory fixed to the indirect maximal center gauge, 
which has the advantage that abelian monopoles can be identified after the first step of 
maximal abelian gauge fixing, and center vortices can be identified after the final center 
projection. 
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We concentrate on static monopoles (monopole current in the time direction) on a given 
time slice. It is found (at jS = 2.4) that almost all monopole and antimonopoles in the time 
slice (97%, on average) lie on P-vortex lines. Along the vortex lines, the predicted alternation 
of monopoles with anti-monopoles is observed in the Monte Carlo data. Moreover, the 
(gauge-invariant) action around a monopole cube in a fixed time slice is highly asymmetric: 
almost all of the excess action is in the plaquettes pierced by a P-vortex! 

■// 

'0^ 



Figure 35: Hypothetical collimation of monopole/antimonopole flux into center vortex tubes 
on the abelian-projected lattice. 
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Figure 36: Excess plaquette action distribution on a monopole cube (left figure) pierced 
by a single P-vortex. For comparison, the excess action distribution is also shown for a 
no-monopole cube (right figure) pierced by a P-vortex. From Ambj0rn et al., ref. [150]. 



The situation found in Monte Carlo simulations (again for the SU(2) group at /9 = 2.4) is 
shown in Figs. 36 and 37. We consider cubes on the lattice, at fixed time, containing static 
monopoles (shaded circles), and which are also pierced by P-vortex lines (thick strip). To 
avoid inhomogeneities due to neighboring monopoles, we also consider "isolated" monopoles 
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Figure 37: Same as Fig. 36, but for "isolated" static monopoles, with no nearest-neighbor 
monopole currents. From Ambj0rn et al., ref. [150]. 



with no nearest-neighbor monopole cubes on the time slice. Define the excess action for 
plaquettes on a monopole cube as 

^ = ^0 - (^TV [C/t/C/^C/1cube plaq) (224) 

where Sq is the usual (gauge-invariant) plaquette expectation value, and the ?7-link variables 
are the hnks on the unprojected lattice. We can calculate S separately for the cases that 
the cube plaquettes are pierced (dark strips), or not pierced (dark lines) by P- vortices 
on the corresponding center projected lattice. What is found is that the excess action is 
overwhelmingly concentrated along the line of the P-vortex (Fig. 36). Moreover, we can 
also calculate the excess action distribution for plaquettes on cubes which are pierced by 
a P-vortex, but which do not contain a monopole. These are also far above the average 
plaquette value. When we look in particular at isolated monopoles (Fig. 37), we find that 
the excess action distribution does not much depend on whether a monopole is, or is not, 
inside a cube. The correlation of excess action is with vortex piercings, not monopoles. 

It is clear that the gauge-invariant field strength around a monopole, far from being 
Coulombic, is instead highly coUimated in the vortex direction. Similar measurements have 
been carried out for cubes which are three or four lattice spacings on a side, and which 
contain either a single static monopole, or no monopole current [150]. Instead of excess 
action, the observable is the fractional deviation of Wilson loop expectation values, for 
loops running around a face of the 3 or 4-cube. The results are qualitatively similar to 
what is seen in Fig. 37 for 1-cubes: There is a large effect for faces of the cube pierced 
by center vortices, and very little correlation with whether or not the cube also contains a 
static monopole. 

Recently the P-plaquette excess-action calculation has been repeated by Gubarev et al. 
[86] at a variety of couplings. These authors confirm that monopole worldlines lie on vortex 
sheets, that the excess action is directed along the vortex, and that vortex density scales. 
Further, they report that the excess action pS is a lattice-independent constant. This 
finding can be interpreted as implying the existence of a sheet of singular action density in 
the middle of center vortices in the continuum limit. 

The strong directionality of field strength along vortex lines whether or not monopoles 
are present, which is found not just on 1-cubes but on larger cubes as well [150], supports 
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the general picture of vortex confinement, in which confining flux is colhmated (at flxed 
time) into tube-hke structures. For isolated monopoles at /3 = 2.4, the monopole locations 
do not appear to play a crucial role as far as the action distribution is concerned. We take 
note of an interesting recent argument by Kovner, Lavelle and McMuUan [155], who point 
out that in SU(2) lattice gauge theory there are residual permutation gauge transformations 
in abelian projection gauges, which can cause monopole-antimonopole pairs to appear or 
disappear along vortex lines. Those pairs, at least, are certainly gauge artifacts. Kovner, 
Lavelle and McMuUan also argue that monopole charge is an ambiguous observable in pure 
SU(2) lattice gauge theory, due to the absence of a charge conjugation operator in the theory. 



8 Coulomb Energy, and the Large N Limit 

In this section I will discuss some aspects of confinement in Coulomb gauge, and in the large 
N limit of SU(N) gauge theory. 



8.1 Confinement in Coulomb Gauge 

The Yang-Mills Hamiltonian in Coulomb gauge has the following form: 

H^ljd'x {E^''^ ■ E^''^ + 5" . 5") + i I d'xd'y roi^)V^\^, y)Joiy) (225) 

where is the zero-component of the conserved color current in eq. (2), is the trans- 
verse color electric field operator, and 

V''\x,y) = \M-\-d'^)M-^Y'' (226) 

L J x,y 

where 

M"^ = -diD^{A) = -d^S"^ - e^^^Aldi (227) 

is the Faddeev-Popov operator, and ^4^ is in Coulomb gauge. Sandwiching the V operator 
between static color sources in representation r, located at points separated by a distance 
i?, defines the static Coulomb potential 

V^^i{R)^{Tv[KV^\Q.R)tl]) (228) 



which is readily seen to be proportional to the quadratic Casimir, since <^V"''y oc 5ab- We 
see that the static Coulomb potential automatically satisfies Casimir scaling. This is not 
surprising, since the static Coulomb potential can also be understood as arising from instan- 
taneous one-gluon exchange in Coulomb gauge. From here on we will consider the Coulomb 
potential for fundamental representation sources in particular. 

There is a scenario for confinement, originally due to Gribov [156], and advocated in 
recent years by Zwanziger [157], which makes essential use of Coulomb gauge and focuses on 
the static Coulomb potential. The idea is roughly as follows: In "minimal" Coulomb gauge, 
the path integral is restricted to configurations such that the Faddeev-Popov operator M 
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has only positive eigenvalues. The boundary of this region in configuration space, where the 
Faddeev-Popov operator acquires negative eigenvalues, is known as the Gribov Horizon. 
Since the dimension of configuration space is very large (on the order of the lattice volume, 
in a lattice formulation), the bulk of configurations will be located close to the horizon (just 
as the volume measure r'^~^dr of a ball in (i-dimensions is sharply peaked near the radius 
of the ball). Because it is the inverse of the M operator which appears in the Coulomb 
potential (228), we may expect that the near-zero eigenvalues of this operator, typical of 
configurations near the horizon, will enhance the magnitude of the Coulomb potential, 
possibly resulting in a confining potential at long distances. 

As Zwanziger has recently pointed out in ref. [158], the Coulomb potential is, at the 
very least, an upper bound on the static quark potential. The argument is rather simple. 
In Coulomb gauge, we can construct a physical state consisting of very massive charged 
sources ("quarks") separated by a distance R 

l*..) = g^"(0)g"(i?)|*o) (229) 

where |\l/o) is the vacuum state and index a denotes the quark color indices. This is, of 
course, not the most general physical state containing two static sources; other states would 
contain also gluon operators acting on the ground state. It does serve, however, as a trial 
state from which one can extract the Coulomb energy. Defining 

Q(i?,t) = gt"(0,t)g"(i?,t) (230) 

the energy expectation value of the state ^^q, above the vacuum energy, is 

[^o\Q^iR,0)HQiR,0)\<iJo] 

= Ese + Vcoul{R) (231) 

where Ege is the quark self-energy, and the i?-dependent part of the energy expectation value, 
for very massive sources, can be identified with the instantaneous Coulomb contribution. 
Although the quark self-energy is divergent in the continuum formulation, it is of course 
regulated on the lattice. 

Next, we note the identity [159] 

S,,^-^rn-^\og[G{R,T)] (232) 

where 

GiR, T) = (q\R, T)Q{R, 0)) (233) 



is a Euclidean vacuum expectation value. In Coulomb gauge, as in any physical gauge, the 
existence of a transfer matrix implies that 

G(it:,r)=^c„exp[-£;Wr] (234) 
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where Ej^^'^ is the energy (above the vacuum energy) of the n-th energy eigenstate having a 
finite overlap with the trial state "ifgg. From this we have the obvious inequality 

= - - (235) 

which is just the statement that the energy expectation value of the trial state is an 
upper bound on the ground state energy of the static quark- antiquark system. This ground 
state energy defines the usual static potential 

= E'^^ + V{R) (236) 

If the static potential V{R) is hnearly rising, then for sufficiently large R the static quark 
self-energies Ege, E'^^ are neghgible compared to V{R), at least on the lattice. Prom this it 
follows that, asymptotically, 

V{R) < V,oui{R) (237) 

Therefore, if the static potential in fundamental representation is confining, the correspond- 
ing Coulomb potential is also confining. Zwanziger and Cucchieri [160] have further sug- 
gested that the bound is saturated, and that the Coulomb potential is the static confining 
potential. If true, this would be a little puzzling, for how would string-like properties 
emerge? If we accept that the confining potential is simply due to instantaneous one-gluon 
exchange (equivalent to the Coulomb energy term in the Hamiltonian) , it is a little hard to 
understand the origin of roughening, or the Liischer term. 

There are a number of different approaches to calculating Vcoui{R) via lattice Monte 
Carlo simulations. One method is to compute the Dqq — (AqAq) component of the gluon 
propagator in Coulomb gauge, using the standard lattice definition 

A,ix) = ^^[U,ix) - Ulix)] (238) 

In the continuum formulation, the instantaneous part of DQo{x,t) is proportional to the 
Coulomb energy. A drawback of this approach is that the modulus of the vector 
potential is bounded in the lattice formulation, and therefore Dqq^R) cannot possibly have 
an asymptotically linear growth on a large lattice.^^ A better method, advocated recently by 
Zwanziger and Cucchieri in [160], is to invert the lattice-regulated Faddeev- Popov operator 
M numerically, in each thermalized lattice configuration, and from that inverse operator 
compute the expectation value on the rhs of eq. (228). The result in momentum space 
{k — |fc|), for k'^Vcoui{k) at /3 = 2.5 in SU(2) gauge theory, is plotted in Pig. 38, together 
with a best fit to 

Bp 

For A ^ 0, the Coulomb potential has linear confining behavior. The intercept of the data 
with the y-axis, which gives the value for A, does indeed appear to be non-zero from the 
data. According to ref. [160], the resulting string tension scales, and appears to nearly 
saturate the inequahty of eq. (237). 

"'^"'^This objection does not apply in Landau gauge, since the gluon propagator in position space is non- 
singular in the infrared [161]. The effect of vortex removal on Landau gauge propagators, and its possible 
implications for confinement, has been investigated recently by Langfeld et al. in ref. [162]. 
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Figure 38: Coulomb energy x|k|^ in momentum space, at f3 = 2.5 in SU(2) lattice gauge 
theory. The solid line is a fit to eq. (239). From Zwanziger and Cucchieri, ref. [160]. 



A third approach gives a somewhat different result, and also indicates a strong relation- 
ship of the Coulomb energy to the center vortex theory. On the lattice, the expectation 
value G{R, T) is simply a correlator of two timelike Wilson lines in Coulomb gauge. Define 



L(x, T) = Uo{x, 0)[/o(x, l)...Uo{x, T - 1) 



so that 



G'(i?,T) = (Tr[Lt(0,T)L(i?,T)]) 
We also define the lattice logarithmic derivative 

V{R,T) =log^ 



(240) 
(241) 

(242) 



G(i?,r + i) 

(with G{R, 0) = 1). It follows that, as the continuum limit is approached at large 

V,oui{R)^V{R,Q)-Ese (243) 

is the Coulomb energy, while 

F(i?)= hm V{R,T)-E',, (244) 

1 — *oo 

is the static quark potential. 



86 



Assuming V{R,T) is asymptotically linear, with string tension cr(T), we must have 



ct(0) Ucoui large (5 

a{T) a large T ^ ' 

Then saturation of the confining potential by the Coulomb potential, a ~ Ocxmii implies that 
(t(T) is approximately T-independent at large /5. This idea has been tested very recently in 
rcf. [163]. In Fig. 39 we display the result for V[R^ 0) at /3 = 2.5, together with two different 
fits to a linear potential (with and without a Liischcr term, whose presence in this case is 
not well motivated). From the figure, the Coulomb potential is clearly linear. However, the 
associated string tension is almost three times that of the asymptotic string tension. In Fig. 
40 we show all of our results for (j(0)/(J, where a is the accepted value for the asymptotic 
string tension. Instead of converging to one, the ratio is actually increasing in the range 
of /5 = 2.2 — 2.5 as (5 increases. This evidence argues that the Coulomb string tension is 
substantially greater than the asymptotic string tension. 

V(R,0)=-log[G(R,1)], |3=2.5 




Figure 39: V{R^ 0) at = 2.5. Open squares are derived from timelike hnk correlators 
1) in Coulomb gauge; open circles are the same data extracted from lattices with 
center vortices removed. From Greensite and Olejnfk, ref. [163]. 

An interesting question, in view of the previous discussion of confinement mechanisms, 
is whether center vortices are also responsible for the Coulomb string tension. This question 
can be addressed by applying the vortex-removal procedure described in section 6.1 above. 
With vortices removed, the modified configuration is fixed to Coulomb gauge by a standard 
over-relaxation method, and T) is computed in the no-vortex ensemble. This result. 
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Figure 40: Ratio of (t(0)/(J, where a{0) —>■ acouiomb in the continuum hmit, and a is the 
usual asymptotic string tension. The string tensions are extracted from fits to a — b/R + aR, 
with either 6 = 0, or the string-motivated value of 6 = 7r/12. Prom Greensite and Olejnik, 
ref. [163]. 



for V{R,0) at /? = 2.5 in the no- vortex ensemble, is also shown in Fig. 39 (open circles). 
Prom this figure it is clear that removing center vortices from the lattice also removes the 
Coulomb string tension completely. 

8.2 Confinement at Large N 

It was pointed out long ago by 't Hooft [164] that the Feynman diagram expansion in SU(N) 
gauge theories could be organized as a double expansion in powers l/g'^N and 1/A^. The 
leading diagrams in for, e.g., a scattering process or a fundamental Wilson loop, are 
planar diagrams; i.e. Feynman diagrams which can be drawn on a plane surface without 
having any two propagators cross one another (except at a vertex). A similar 1/N organi- 
zation exists in the lattice strong-couphng expansion. It is hoped that the case oi N — 2 or 
N = 3 can be regarded as a relatively small perturbation of the N = oo limit. Perhaps this 
limit is more tractable analytically; perhaps it leads us to new insights. 

There is one striking property of the N = oo limit which seems quite unlike the situation 
at N — 2, 3; this is the property of factorization. Let Oi and O2 be any two gauge-invariant 
operators. Then to leading order in 

{0,02) = (Oi) (O2) (246) 

Factorization, applied to Wilson loops, implies that Casimir scaling is exact in the — > 00 
limit, as already noted in section 4.2 above. But this property also implies, since AO^ = 

^^One might say that in this hmit, the Casimir scahng region in Fig. 31 expands, and the N-ahty dependent 
region recedes, out to infinite distances. 
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(02) - {0)\ that the rms deviation of any gauge-invariant operator from its mean value 
vanishes as — > oo. This has an astonishing consequence. If we would imagine performing 
a numerical simulation of SU(N) lattice gauge theory at some enormous value of N , such 
that all non-leading powers of N could be neglected, then every thermalized configuration 
would give almost the same value for any given observable. This led Witten [40] to propose 
the idea of a large-N master field; i.e. a single configuration in the N = oo limit 

with the property that for any gauge-invariant functional 0[A] of the gauge field, 

(O) = 0[A;^"^^'=''] (247) 

The master field is certainly not unique, since any thermalized configuration will do. There 
exists, in fact, a fairly simple equation whose solution is known to yield QCD master field 
configurations [165]. Unfortunately, although the equation can be solved perturbativcly and 
reproduces the planar Feynman diagram results, it seems to be no easier to solve at the 
non-perturbative level than any other formulation of QCD. 

SU(2) and SU(3) lattice gauge theory certainly do not come close to having the AO = 
property for Wilson loops. Consider the frequency distribution of Tr [[/(O)], which can be 
obtained in Monte Carlo simulations by plotting a histogram for the values obtained for 
Tr[t/(0)] around loop C in thermalized configurations. In SU(2) lattice gauge theory, for 
loops whose extension is comparable to the confinement scale, the probability distribution 
is very fiat, as sketched in Fig. 41. The exponentially small value of the Wilson loop 
is obtained in the average not because the loop holonomy typically has a small trace, but 
rather because positive and negative values of the trace cancel almost entirely in the average 
over configurations. In contrast, the frequency distribution at N = oo would have to be a 
delta function. We must clearly be cautious in extrapolating N = oo physics to A/" = 2 or 

Nevertheless, in some ways SU(N) gauge theory at small N does seem to be close to the 
N = oo limit, in the sense that ratios of physical quantities converge rather rapidly as N 
increases. This fact can be seen in Fig. 42, which shows the results of Lucini and Teper 
[166] for the masses of the lowest-lying 2++, and first excited 0++ glueballs, divided 
the square root of the string tension, for = 2 — 5. Figure 43 displays the string tension 
as a function of the tadpole-improved lattice version of the 't Hooft coupling g'^N 

where /3 is the usual Wilson action lattice couphng. At fixed lattice 't Hooft couphng, there 
is obviously not much variation in the string tension for different A", indicating that there 
are only small deviations from the planar limit, at least for this quantity. 



8.2.1 k-string Tensions at Medium N 

k-string tensions are the asymptotic string tensions of the lowest dimensional represen- 
tation with N-ality k, represented by a single column of k boxes in the Young tableau. 
Because of color screening, a /c-string tension is also the asymptotic string tension of any 
color representation of N-ahty k. 
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Figure 41: Schematic Wilson loop frequency distributions for a) SU(2), and b) SU(oo). The 
SU(2) distribution is peaked at the center elements, and the finite value of Wilson loops 
is due to a tiny asymmetry between Tr (^f) and Tr(— ^f). At iV = oo, every evaluation of 
ReTr[t/(C)] yields the same number. 
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Figure 42: Masses of the lowest lying scalar, tensor, and first excited scalar glueballs, from 
lattice simulations with = 2 — 5 colors, in units of i/a. The m/y/a ratio is plotted vs 
l/N"^. From Lucini and Teper, ref. [166]. 



Douglas and Shenker [44], generalizing the work of Seiberg and Witten [152] to the SU(N) 
gauge symmetry, derived from softly broken TV = 2 SUSY a supersymmetric version of the 
dual abelian Higgs model, in which the (dual) U{1)^~^ symmetry^roup is spontaneously 
broken. In the dual theory there are distinct chiral superfields M„, Mn, n = 1, ..,N — 1 with 
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Figure 43: Square root of the lattice string tension vs. the (tadpole-improved) lattice 't Hooft 
coupling A/ = gjN. From Lucini and Teper, ref. [166]. 



magnetic monopole components, and vacua at 

{MM)n oc sin[7m/iV] (249) 

This in turn leads to iV — 1 distinct string tensions T„ associated with each of the N — 1 
magnetic U{1) symmetries, with 

T„ oc sm[irn/N] (250) 

Let q"" be the color components of a quark field in the fundamental representation of SU(N). 
The product of components q^q^q^ ...q^ , which has N-ality k, is charged only under the k-th 
magnetic U(l) factor, and is neutral with respect to the other factors [45]. The upshot is 
that the representation of lowest dimension of N-ality k, which is an antisymmetric product 
of k q'^ factors, would have an asymptotic string tension (y{k) = T^. This is the Sine Law 
prediction for the N-ality dependence of asymptotic string tensions 

^^^'^) = ^ = ^hIRiVr (251) 

The same dependence is predicted in a related M-theory version of supersymmetric gauge 
theory known as MQCD [45]. 

For comparison, suppose that Casimir scaling were exact in the intermediate distance 
regime. The ^-representations cannot be screened by gluons, so it may be that for these 
representations the string tension in the Casimir scaling regime is also the string tension in 
the asymptotic regime. Under those two assumptions, the Casimir scaling prediction is that 

R(k,N)^'<§^ P52) 

With this motivation, the authors of refs. [56] and [167] have calculated R{k^ N) for the 
SU(4) and SU(6) gauge groups, for ordinary non-supersymmetric lattice gauge theory, with 
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the results shown in Table 1. For D = 2 + 1 dimensions, Casimir scaling seems preferred, 
while in D = 3 + 1 dimensions there appears to be remarkable agreement with Sine Law 
scahng.^^ 



Table 1: Casimir Scaling and Sine Law Ratios vs. data in 2+1 and 3+1 dimensions. 



R(k,N) 


Cas Sin 


Data, D = 2 + 1 [56] 


Data, D = 3 + 1 [167] 


R(2,4) 
R(2,6) 
R(3,6) 


1.33 1.41 
1.60 1.73 
1.80 2.00 


1.355(6) 
1.616(9) 
1.808(25) 


1.403(15) 
1.72(3) 
1.99(7) 



The interpretation of these results is a little obscure, in part because the Sine scaling law 
is non-universal, and is subject to 1/iV^ corrections even in the case of softly broken Af = 2 
SUSY [168]. There is no obvious reason that the Sine Law should hold very accurately 
at, say, N=4, especially in a non-supersymmetric theory. Nor does the supersymmetric 
case cast much light on the approximate Casimir scaling of the non-supersymmetric theory, 
where the string tension at intermediate distances depends on the quadratic Casimir of the 
color representation, rather than the abelian electric quark charges. 

In terms of the vortex picture, if we assume short-range correlations of center vortex 
fields in a plane (eq. (94)), and assume cither Casimir or Sine-Law scaling for the k-string 
tensions, then it is not hard to work out the vortex densities at a given [62]. These have 
no obvious pathologies in the large-N limit, and in particular do not grow with N. 

As N increases, the Casimir and Sine scaling laws converge to a common limit for the 
k-string tensions, ior N k 

R{k, N)^k (253) 

It is interesting that this limit has a simple explanation in the context of the center vortex 
mechanism [62] . Stated briefly, the confining dynamics at asymptotic distances are described 
at large iV by a center monopolc Coulomb gas, in which the analogue of electric charge is 
N-ality. We recall that for any N > 2 there exist center monopoles as well as center vortices; 
these monopoles are the joining point of N A; = 1 vortices, and go over to U(l) monopoles in 
the N ^ oo limit. The /c-string tension in a monopole Coulomb gas is proportional to the 
electric charge, which in this case is k, because there are k independent fiux tubes stretching 
between the quark and antiquark. 

8.2.2 The QCD String at Large N 

One of the reasons for being interested in the large N limit is that this limit suggests a string 
description of the QCD flux tube. The possible interpretation of high-order planar diagrams 

^•^Thc results for D = 3 + 1 dimensions for the SU(5) gauge group, reported in ref. [56], are not in such 
good agreement with the Sine Law, but the results of ref. [167] are claimed to be more accurate, so that is 
what we quote here. 
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as string worldsheets, as indicated in Fig. 44, was pointed out by 't Hooft in his seminal 
article on large N gauge theory [164]. In more recent years, the AdS/CFT conjecture [3] 
has actually provided a precise string representation of Wilson loop expectation values in 
the N — oo limit, at least for large 't Hooft couphngs g'^N » 1. On the other hand, 
the string of the AdS/CFT approach is not quite the same thing as the QCD flux tube. 
For one thing, the AdS/CFT string lives in 10 dimensions, and for another, there is also 
an AdS/CFT string representation for Wilson loops in the non-confining case of unbroken 
Af — A supersymmetric gauge theory, where there is no flux tube at all. 

Returning to the planar diagram- as- worldsheet idea, note that a time-slice of a high- 
order planar diagram for a Wilson loop reveals a sequence of gluons, with each interacting 
only with its nearest neighbors in the diagram (Fig. 44). This picture suggests that the 
QCD string might be regarded, in some gauge, as a "chain" of gluons, with each gluon 
held in place by its attraction to its nearest neighbors in the chain. The linear potential 
in this gluon chain model comes about in the following way: As heavy quarks separate, 
we expect that at some point the interaction action energy increases rapidly due to the 
running coupling. Eventually, it becomes energetically favorable to insert a gluon between 
the quarks, to reduce the effective color charge separation. This fits in nicely with the 
result of ref. [163], where it is found that the Coulomb force of a quark-antiquark state, 
with no constituent gluons, is much higher than the accepted asymptotic force. As the 
quarks continue to separate, the process repeats, and we end up with a chain of gluons. 
The average gluon separation R along the axis joining the quarks is fixed, regardless of the 
quark separation L, and the total energy of the chain is just the energy per gluon times the 
number of gluons in the chain, i.e. 

E chain ~ IT' gluons E gluon — — ^ — L — gL (254) 

In this picture, the linear growth in the number of gluons in the chain is at the heart of the 
long-range static potential. 

The gluon-chain picture, originally proposed many years ago by Thorn [169] and the 
author [170, 159], and recently revived in rcf. [171], has quite a number of nice features. 
For one thing, because of its large-N origin, the gluon chain model naturally accounts 
for Casimir scaling in the large N limit. For example, the string tension in the adjoint 
representation is twice that of the fundamental representation simply because there are 
two gluon chains stretching between sources in the adjoint representation. But also the 
correct N-ality dependence is obtained by (l/N"^ suppressed) string-breaking interactions. 
As discussed in detail in rcf. [171], the model also naturally accounts for the logarithmic 
growth of the QCD flux tube with quark separation (roughening), and the existence of a 
Liischer term. 

The gluon chain model is a representation of the QCD flux tube in terms of particle (i.e. 
gluon) excitations, and is in some sense "dual" to the description of the QCD vacuum in 
terms of field fluctuations (e.g. vortices and monopoles). It is hoped that this representa- 
tion of the QCD string, in which gluon separations (and therefore effective couplings) are 
bounded, may lend itself to variational treatments of the QCD string and low-lying glueball 
states, as outhned in ref. [171]. Another possible direction for analytic work is the approach 
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Figure 44: The gluon chain as a time shce of a planar diagram (shown here in double-line 
notation). The solid (open) hemisphere denotes a quark (antiquark) color index. 



initiated by Bardakci and Thorn [172], in which planar diagrams in light-cone gauge can 
be explicitly reformulated as string amplitudes. In Coulomb gauge, the gluon-chain picture 
lends itself to tests via Monte Carlo simulations, some of which were reported long ago in 
the second article of ref. [170]. 



9 Conclusions 

I have emphasized center symmetry in this article because there are good reasons to believe 
that center symmetry is important. To review briefly: The confinement and deconfine- 
ment phases at finite temperature are the symmetric and the spontaneously broken phases, 
respectively, of a global Zn center symmetry. The existence of a finite string tension is 
related, by rigorous arguments, to the behavior of the center vortex free energy. Even in 
theories with a trivial center subgroup, such as SO (3) lattice gauge theory, the vortex free 
energy is a good order parameter. When global Zj^ center symmetry is explicitly broken 
by scalar fields, the phase transition from a Higgs to a distinct confinement phase is lost. 
Separate Higgs and confinement phases can exist only when the scalar fields are in zero 
N-ality representations, preserving the global center symmetry of the Lagrangian. Finally, 
the asymptotic string tension of static quarks depends on their color charge only through 
the transformation properties of the quarks under the center subgroup. 

In view of these facts, a confinement mechanism based on center vortices would seem 
to be the natural way of realizing global center symmetry in the vacuum state. In recent 
years it has become possible to identify vortex locations in lattices generated by Monte 
Carlo simulations, via center projection in an adjoint gauge. Identification of objects in the 
vacuum by a gauge-fixing procedure is always subject to ambiguities, but these worries can 
be at least ameliorated by demonstrating the correlation of vortices with gauge- invariant 
observables. In particular: The number of vortices mod 2, linked to a large Wilson loop in 
SU(2) lattice gauge theory, is found to be strongly correlated with the sign of the Wilson 
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loop. The vortex surface is associated with an average (gauge-invariant) plaquette action 
which is well in excess of the vacuum average. Removing vortices from unprojected lattice 
configurations removes the string tension, removes chiral symmetry breaking, and sends 
every configuration to the trivial topological sector. Apart from these correlations with 
gauge invariant observables, it is found that vortices closed by lattice periodicity in time 
explain the existence of a spatial string tension across the deconfinement transition. The 
density of center vortices scales according to asymptotic freedom, and the string tension of 
the projected lattice accounts, pretty nearly, for the accepted asymptotic string tension. 

The abelian monopole theory, based on abelian projection, can also claim some numerical 
success. Monte Carlo investigations of vortices and monopoles should not be irreconcilable; 
I believe that the relevant numerical fact is that abelian monopoles lie along vortex lines 
in a monopole-antimonopole chain, as discussed in section 7. In general, abelian monopole 
mechanisms are most attractive in models which contain an adjoint representation scalar 
field in the Lagrangian, and therefore naturally single out a unique abelian [/(l)^"^ sub- 
group. QCD, on the other hand, does not single out a unique abclian-projcction subgroup, 
and the Casimir scaling of string tensions found at intermediate distances, together with 
N-ality dependence asymptotically, argue against the idea of an underlying U{1)^~^ mech- 
anism. It is possible, however, that monopole worldlines on vortex sheets do play a role in 
the generation of topological charge. 

We have some pieces of the puzzle, but they are not all in place. It seems very likely 
that center symmetry and center vortices are an important part of the picture. But despite 
the numerical evidence obtained in recent years, and despite the input from supersymmetry 
and M-theory, quark confinement in QCD is imperfectly understood, and an actual proof is 
elusive. The confinement problem is still open, and remains a major intellectual challenge 
in our field. 
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